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ABSTRACT 

We present a complete sample of 24 radio-loud quasars (RLQs) from the new 
7C Redshift Survey. Every quasar with a low-frequency (151 MHz) radio flux-density 
5i5i > 0.5 Jy in two regions of the sky covering 0.013 sr is included; 23 of these have 
sufficient extended flux to meet the selection criteria, 18 of these have steep radio 
spectra (hereafter denoted as SSQs). The key advantage of this sample over most 
samples of RLQs is the lack of an optical magnitude limit. By combining the 7C and 
3CRR samples, we have investigated the properties of RLQs as a function of redshift 
z and radio luminosity £151. 

We derive the radio luminosity function (RLF) of RLQs and find that the data are 
well fitted by a single power-law with slope a± = 1.9±0.1 (for H — 50 km s _1 Mpc _1 , 
fiM = 1, Oa = 0). We find that there must be a break in the RLQ RLF at log 10 (Li5i/ 
W Hz -1 sr -1 ) < 27, in order for the models to be consistent with the 7C and 6C source 
counts. The z-dependence of the RLF follows a one-tailed gaussian which peaks at 
z = 1.7 ± 0.2. We find no evidence for a decline in the co-moving space density of 
RLQs at higher redshifts. 

A positive correlation between the radio and optical luminosities of SSQs is ob- 
served, confirming a result of Serjeant et al. (1998). We are able to rule out this 
correlation being due to selection effects or biases in our combined sample. The radio- 
optical correlation and best-fit model RLF enable us to estimate the distribution of 
optical magnitudes of quasars in samples selected at low radio frequencies. We con- 
clude that for samples with 6*151 < 1 Jy one must use optical data significantly deeper 
than the POSS-I limit (R w 20), in order to avoid severe incompleteness. 

Key words: Galaxies: active - quasars: general - quasars: emission lines 



1 INTRODUCTION 



The variation of the space densities and properties of quasars 
with redshift contains information about structure forma- 
tion and evolution at different epochs. Radio galaxies are 
important for similar reasons and their relationship with 
radio-loud quasars (RLQs) is now under scrutiny. Accord- 
ing to orientation-based unified schemes, RLQs and powerful 
radio galaxies are the same type of objects, their apparently 
different properties being due to a systematic difference in 
the angle between our line-of-sight and the radio jet axis (e.g. 
Scheuer 1987, Antonucci 1993). In orientation-based unifi- 
cation, the broad line region and quasar nucleus are hidden 
in radio galaxies behind an optically-thick torus. If there is 
not a sharp cut-off in the optical depth of the torus, then 
objects close to the transition angle between radio galaxies 



and quasars may appear to be faint reddened quasars (Baker 
1997). 

One of the major problems in trying to understand the 
collective properties of a sample of objects is ensuring that 
the sample provides reliable information on the parent pop- 
ulation. Particularly in astrophysics, care must be taken to 
minimise and understand selection effects. For example, be- 
cause of the steepness of the radio luminosity function, there 
is a strong correlation between radio luminosity and redshift 
in the 3CRR sample (Laing, Riley & Longair 1983, here- 
after LRL). Hence, it is impossible to disentangle correla- 
tions with radio luminosity from those with evolution with a 
single flux-limited sample. Radio quasar samples selected at 
high-frequency contain mostly flat-spectrum quasars (FSQs, 
«rad < 0.5 Ij), where a ra( j is evaluated at 1 GHz in the 
source rest-frame. The radio emission from these objects is 
often dominated by optically-thick radio cores, which are en- 
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hanced by Doppler boosting because their jet axes lie close 
to the line-of-sight. Therefore their total radio flux-densities 
(and hence probability of being above a given flux-limit) are 
orientation-dependent and often strongly time variable (e.g. 
Seielstad, Pearson & Readhead 1983). In contrast, the radio 
emission from steep-spectrum quasars (SSQs, a rad > 0.5) 
is dominated by optically-thin extended emission and hence 
low-frequency selected samples suffer far less greatly from 
these orientation-dependent or variability biases. Some radio 
sources, particularly compact-steep spectrum (CSS) sources 
(e.g. Fanti et al. 1995) and Giga-Hertz Peaked Spectrum 
(GPS) sources (e.g. O'Dea, Baum & Stanghellini 1991), suf- 
fer synchrotron self-absorption at low-frequencies and hence 
have radio spectra which flatten at low-frequency, eg. 151 
MHz. By evaluating the spectral indices at rest- frame 1 GHz 
we avoid rejecting these source from our sample, although 
the sensitivity of a survey to such sources is clearly a function 
of both survey frequency and source redshift (see Blundell, 
Rawlings & Willott 1998a (Paper IV) for a fuller discussion). 

The luminosity function describes the co-moving den- 
sity of a class of objects as a function of redshift and lumi- 
nosity. For quasars, one can define the luminosity function 
in terms of the optical luminosity function (OLF), X-ray lu- 
minosity function (XLF), or the radio luminosity function 
(RLF). Most large quasar surveys have been selected us- 
ing optical, and more recently X-ray, criteria so the quasar 
OLF and XLF are better constrained than the RLF. Boyle, 
Shanks & Peterson (1988) showed that the OLF can be fitted 
with a broken power-law with the break luminosity evolving 
with redshift — an effect normally termed pure luminosity 
evolution (PLE). However, Hewett, Foltz & Chaffee (1993) 
and Miller et al. (1993) showed that PLE was an inadequate 
model to explain the evolution of luminous quasars. Hawkins 
& Veron (1995) explained the break in the power-law of 
Boyle et al. as due to incompleteness in their data. Hewett 
et al. proposed that the power-law index of the OLF steep- 
ens as the redshift increases. This is supported by the recent 
analysis of the quasar OLF by Goldschmidt & Miller (1998). 
The space density of the most luminous quasars shows lit- 
tle or no evolution (Miller et al. 1993), whilst the optically 
fainter quasars show significant evolution (a factor of ~ 50 
increase between z ~ 0.5 and z ~ 2). The XLF also appears 
to have a broken power-law form (Maccacaro et al. 1991; 
Delia Ceca et al. 1992). Boyle et al. (1994) find the best- 
fitting model for the XLF is pure luminosity evolution for 
z < 1.6 and a constant co-moving density at higher redshift. 
However, it must be remembered that these quasar samples 
selected in the optical and X-ray wavebands may be prone 
to selection effects. If, for example, many quasars are highly 
reddened (Webster et al. 1995), then this would greatly alter 
the shape and evolution of the OLF, and because of absorp- 
tion by neutral hydrogen, the XLF. 

One of the difficulties of determining the quasar RLF is 
that radio surveys typically contain mostly radio galaxies, 
so all sources must be reliably identified as quasar or radio 
galaxy. Dunlop & Peacock (1990) derived the RLF sepa- 
rately for flat- and steep-spectrum radio sources. However, 
the high selection frequency (2.7 GHz) of the samples used 



* The convention for spectral index, ct ra d, is that S v ctv c "rad j 
where S v is the flux density at frequency v. 



in this work meant a large fraction of flat-spectrum sources 
in their study and hence the steep-spectrum RLF was more 
poorly constrained. The steep-spectrum RLF fit of Dunlop 
& Peacock includes both radio galaxies and steep-spectrum 
quasars. No attempt was made to separate these different 
types of object. Hence the steep-spectrum quasar RLF has 
not previously been well-constrained. Radio selected quasar 
samples usually have additional optical selection criteria, 
which can be a cause of incompleteness and selection biases. 
Since the 7C quasar sample has no optical selection criteria 
it is complete, so in conjunction with the 3CRR sample of 
LRL it can be used to derive the RLF for the low-frequency 
selected population which is dominated by SSQs. 

There have been several suggestions of correlations 
between the radio and optical properties of AGN. The 
close relationship between narrow-line luminosity and low- 
frequency radio luminosity (e.g. Baum & Heckman 1989; 
Rawlings & Saunders 1991; Willott et al. 1998b (Paper III)) 
suggests a physical link between these two properties. Cor- 
relations between the radio and optical luminosities of RLQs 
have been found (e.g. Browne & Murphy 1987) but incom- 
pleteness worries mean that there are unknown selection ef- 
fects which could affect the results. Correlations between the 
radio and optical may be due to orientation-dependent emis- 
sion (e.g. Baker 1997). This is a problem for flat-spectrum 
RLQ samples because often the radio luminosities are domi- 
nated by a Doppler-boosted core and the optical luminosities 
too may therefore be orientation-dependent (e.g. Wills & 
Lynds, 1978; Browne & Wright 1985). Serjeant et al. (1998) 
presented evidence for a radio-optical correlation for SSQ 
samples. They used the 3CRR sample and a sample based 
on cross-matching the low-frequency (408 MHz) MRC radio 
catalogue with APM data - the MAQS sample. This sample 
has an optical magnitude limit as well as a radio flux limit. 
They also calculated B-band absolute magnitudes using ob- 
served B magnitudes and if-corrections. For high-redshift 
quasars, this radiation was emitted in the UV and so will 
be greatly diminished if there is significant reddening of the 
quasar spectrum by dust. 

Excepting 3CRR, all other SSQ samples also have op- 
tical magnitude limits, usually in the blue part of the spec- 
trum (e.g. Vigotti et al. 1997; Riley et al. 1998). Hence faint 
(and/or reddened) quasars can be missed in these samples. 
If a substantial fraction of quasars are reddened then quasar 
samples with bright magnitude limits will miss them and an 
unrepresentative sample of the quasar population (i.e. only 
the 'naked' optically-bright ones) is seen. The only way to 
be sure of avoiding these selection effects is to identify all the 
radio sources in a flux-limited sample. However, this can be 
very time consuming because most radio samples are domi- 
nated by radio galaxies, not quasars. 

Here we present a complete sample of 24 quasars from 
the new 7C Redshift Survey. In this survey we have obtained 
optical spectra of every radio source in two patches of sky 
(whose total area covers 0.013 sr), with a 151 MHz flux den- 
sity £151 > 0.5 Jy. We made high-resolution VLA maps of 
these sources to ensure reliable optical identifications. All 77 
radio sources have optical / near-infrared identifications and 
only 6 lack secure spectroscopic redshifts. All 6 of these are 
galaxies with estimated redshifts from multi-colour photom- 
etry. Hence we can be confident that there is no significant 
incompleteness in our sample of radio-loud quasars. 
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The layout of this paper (Paper II) is as follows. In 
Section 2 we present the data on the complete 7C quasar 
sample. In Section 3 we discuss our definition of a quasar 
in order to ensure that we are comparing similar types of 
objects from the 3CRR and 7C samples. This is necessary 
because the traditional definition of a quasar is no longer 
tenable in the light of unified schemes. The RLQ radio lu- 
minosity function is derived in Section 4 and compared with 
the optically-selected quasar OLF. In Section 5 we discuss 
the evidence for a correlation between the optical and radio 
luminosities of steep-spectrum quasars. In Section 6 we use 
our knowledge of the radio-optical correlation and the RLQ 
RLF to estimate the completeness of quasar samples selected 
with low- frequency radio and optical flux limits. The opti- 
cal properties of the 7C quasar sample are to be discussed 
elsewhere (Willott et al. in prep., Paper VIII). Also in this 
forthcoming paper, we investigate the amount of reddening 
in the 7C quasar sample and its impact on the completeness 
of quasar surveys. 

We assume throughout that H a = 50 km s _1 Mpc -1 . All 
the calculations in this paper were done for three different 
cosmological models to account for the current uncertainty 
in the values of the matter density parameter (Qm) and the 
cosmological constant parameter (Oa). The models consid- 
ered are I: Q M = 1, fi A = 0; II: Q M = 0, f2 A = 0; III: 
Q M = 0.1, Ov = 0.9. 



2 THE 7C REDSHIFT SURVEY 

Only brief details of the 7C Redshift Survey are given here; 
full details will appear elsewhere (see also Rawlings et al. 
1998; Willott et al. 1998a; Blundell et al. 1998b, 1998c 
(Paper I). The 7C survey was carried out with the Cam- 
bridge Low Frequency Synthesis Telescope (CLFST) at 151 
MHz (McGilchrist et al. 1990) at a resolution of 70 x 70 
cosec(Dec) arcsec. The 7C Redshift Survey was pursued in 
three separate patches of sky: 7C-I, 7C-II and 7C-III. In 
this paper we consider the 7C-I and 7C-II regions, which 
were chosen by Rossitter (1987) to overlap with the 5C6 and 
5C7 fields of Pearson & Kus (1978), respectively. The 7C- 
III region is centred on the North Ecliptic Pole (overlapping 
with the 38 MHz selected sample of Lacy et al. 1993) and 
has similar spectroscopic completeness to 7C-I and 7C-II. 
However the 7C-III objects lie above the declination limit 
of the UKIRT and we have been unable to obtain K-band 
images of most of these objects. For this reason the 7C-III 
region will not be considered further in this paper. The 7C- 
I and 7C-II 151 MHz data were subsequently re-analysed 
(and in the case of the 5C6 field re-observed) by Julia Riley. 
VLA data, where available, were retrieved from the archive 
(by kind permission of Barry Clark), or were recently ob- 
tained by us: full details of all these radio observations and 
the sample selection are described by in Paper I. The 7C-I 
region covers an area of 0.0061 sr and contains 37 sources 
with 5*151 > 0.51 Jy. The 7C-II region covers an area of 
0.0069 sr and contains 40 sources with 5*151 > 0.48 Jy. All 
the sources in the final sample have now been mapped with 
the VLA to produce high-resolution radio maps. These are 
necessary to resolve the lobes and core and get an accurate 
position for the expected optical counterpart. The angular 
separations between the 7C radio source positions and the 



optical positions of the quasars are listed in Table 1. Note 
that these separations are less than 10 arcsec in all cases 
(see Riley et al. 1998 for further discussion of this point). 

Despite the low selection frequency of the 7C survey 
(151 MHz), it is possible that there could be Doppler- 
boosted quasars in the sample whose extended flux alone 
would not exceed the sample flux limit. By careful spectral 
fitting of the data, using integrated flux-densities at a num- 
ber of frequencies from many radio surveys and our own 
VLA maps, we evaluated 1 GHz rest-frame spectral indices 
(hereafter a ra d) for all the sources (as described in Blundell 
et al. 1998a). By consideration of these radio spectra and 
by examining our own VLA maps (Paper I), we were able 
to identify any 7C sources which exceeded the flux limit 
only because of the contribution of their Doppler boosted 
cores. Although six of the 24 quasars in the 7C sample have 
«rad < 0.5, we believe that just one (5C7.230) of these is 
promoted into the sample by such orientation-dependent ra- 
dio emission. The other five quasars in the 7C sample with 
«rad < 0.5 we believe belong to a class of objects we term 
Core- JetS sources (CJSs); the details of their radio structure 
and spectra are discussed in Paper I, but a summary of their 
properties is as follows: these are quite compact sources, 
whose spectra may be showing evidence of synchrotron self- 
absorption, which when resolved with MERLIN or VLBI 
reveal jet-type structures typically on both sides of the core. 

We will use 23 7C RLQs in much of the subsequent 
analysis, comprising 18 which are SSQs (with ow > 0.5), 
and 5 which have extended flux which exceeds the flux-limit 
but which have a ra d ~ 0.4. 

A similar analysis of the 3CRR sample (see Appendix 
B) leads to a sample of RLQs consisting of 39 SSQs and 
one CJS (3C286). Possible reasons for the differing SSQ:CJS 
ratios in the two samples will be explored elsewhere. Two 
quasars from 3CRR (3C345 and 3C454.3) are excluded be- 
cause they are promoted into the sample by Doppler boosted 
core emission. 



2.1 Source identification and if-band imaging 

All 7C-I and 7C-II objects have been imaged at K-b&nd 
(2.2pim) for as long as required to obtain good signal-to- 
noise on the radio source identification (except the Doppler- 
boosted flat-spectrum quasar, 5C7.230). These images and 
their analysis will be published elsewhere (Willott et al, in 
prep), but a log of the observations is given in Table 1. Find- 
ing charts for the entire 7C Redshift Survey were generated 
from the APM catalogue at the RGO, Cambridge (Irwin, 
Maddox & McMahon 1994). This catalogue utilises scanned 
plates from the Palomar Observatory Sky Survey (POSS). 
The astrometry of all images was fixed with the IRAF gasp 
package using these finding charts. K magnitudes were mea- 
sured with the IRAF apphot package using circular aper- 
tures of 5 arcsec diameter where possible. The infrared coun- 
terpart of the radio source was observed in every single im- 
age. Hence, we are confident that we have completely iden- 
tified all members of the 7C Redshift Survey. The measured 
quasar if-magnitudes are listed in Table 3. 
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Name 


Date 


Exposure 
lime (s) 


Seeing 
(arcsec) 


Optical/NIR Position 
(B1950.0) 


7C— optical 
separation (arcsec) 


5C6.5 


971an,30 


240 


1.0 


02 04 14.40 


+33 22 20.1 


3.4 


5C6.8 


971an30 


240 


1.0 


02 04 39 63 


+31 37 51 7 

TOl Kit *J -L • f 


1.6 


'SCfi 33 


971an29 


540 


1.0 


02 08 36 32 


4-33 34 18 5 

| -J • } • >i 1 » ' . 7 


2.8 


5C6.34 


96Scr>04 


540 


1.0 


02 08 38 50 


4-31 21 30 2 

| - ) 1 — 1 • > V7 . — 


3.6 


KCf] 1Q 


9fi lan21 


1080 


1.1 


02 08 53 76 

WO tjij. I u 


4-S2 14 46 1 


2.5 


5C6 95 


96Ser>04 


540 


1.1 


02 11 49.24 


4-29 36 40 


8.3 


KCf] IfiO 


971an29 


540 


0.9 


02 1 4 44 80 


4-33 21 1 5 


0.9 


5C6 237 


971ari29 


540 


1.1 


02 17 49 86 


4-32 27 23 3 


2.0 


5C6 251 


971an29 


540 


1.1 


02 1 8 50 44 


4-34 02 38 9 


4.7 


5C6 264 


971an29 


540 


0.9 


02 19 49 07 


4-S3 43 41 3 


2.2 


'iCfi 282 


96Ser>04 


540 


1.0 


02 21 40.78 


4-34 06 09 2 

| ■ ) i V 7 V7 V ' -7 . — 


0.7 


5C6 286 


971an30 


360 


1.1 


02 22 09.32 


4-S1 45 49 2 


4.4 


5C6 287 


96Ser>04 


540 


0.9 


02 22 19.95 


4-31 39 00 6 

|-7 1 ' / • 7 \ ' V7 . V 7 


4.5 


5C6 288 


95Mar01 


1080 


1.1 


02 22 19.33 


4-31 05 27 5 


1.1 


5C6 2Q1 


971an29 


1080 


1.1 


02 23 09 74 


4-34 08 00 6 


0.6 


7008084-2854 

1 v_ w 1 1 V7 1. ' (7 ' 7 i 


971an30 


540 


1.2 


08 08 32 11 

V / 17 V ) l_J -.7 . 1 1 


4-28 54 01 8 


5.7 


5C7.17 


£7V71V±<Jj! J_l_7 


540 


1.1 


08 09 56.36 


4-27 00 57.5 


6.2 


5C7.70 


971an30 


540 


1.3 


08 14 00.93 


4-29 27 30.9 


1.3 


5C7.85 




540 


1.2 


08 14 45.58 


4-27 26 00.7 


1.1 


5C7.87 


96MarlO 


540 


1.2 


08 15 03.89 


4-24 58 27.3 


6.0 


5C7.95 


95Mar01 


1080 


1.3 


08 15 20.01 


4-24 45 05.4 


6.8 


5C7.118 


96MarlO 


540 


1.2 


08 16 15.06 


4-26 51 29.3 


1.7 


5C7.194 


96MarlO 


540 


1.0 


08 19 14.39 


4-25 48 09.8 


0.9 


5C7.195 


95Mar01 


2520 


1.0 


08 19 20.11 


4-25 15 37.8 


0.8 


5C7.230 








08 21 34.26 


4-24 48 28.5 


3.0 


7C0825+2930 


97Jan30 


240 


1.3 


08 25 05.45 


4-29 30 17.9 


3.2 



Table 1. Log of if-band observations of the 7C quasars and broad-line radio galaxies using IRCAM3 on the UKIRT. The observations 
on 96Sep04 and 96Marl0 were part of the UKIRT service observing program. All other observations were made by us. The 7C-optical 
separation is the angular separation (in arcsec) between the optical position and the radio position given in the 7C catalogue. The 
position of 5C7.230 is from its identification in the APM catalogue, since no A"-band imaging of this object has been obtained. 



2.2 Optical spectroscopy 

About half of the 7C-I/7C-II quasars (13 out of 24) were 
identified on the POSS-I plates as the likely optical coun- 
terpart of the radio source, prior to spectroscopy. Some of 
the other quasars were identified from prior imaging as can- 
didate identifications. Where the object was identified, an 
optical spectrum was taken at the position of this object. If 
an object was not identified before time with a spectrograph 
was available, then a spectrum was taken using blind offset- 
ting to the radio position, aligning the slit with the radio 
axis (Rawlings, Eales & Warren 1990). Most spectra were 
obtained with the ISIS long-slit spectrograph at the WHT. 
Table 2 lists the details of the quasar spectroscopy. The at- 
mospheric conditions in January 1994 were not photometric 
since clouds affected the transparency. 

The spectra were reduced in the IRAF package by sub- 
tracting the bias, flat-fielding with a combination of twi- 
light sky and tungsten lamp spectra, wavelength calibra- 
tion with comparison arc lamps, extinction correction, back- 
ground subtraction and flux-calibration using spectrophoto- 
metric standards. One-dimensional spectra were extracted 
from the 2-D frames using a FWHM aperture for a good 
signal-to-noise ratio (snr) and a FWZI aperture for flux 
measurements. Atmospheric absorption was corrected for by 
scaling the absorption from a standard star spectrum by the 
airmass of the quasar observation. Distortions where the red 



and blue spectra meet (due to the dichroic) were also cor- 
rected using the flux-standard spectra. Final wavelength cal- 
ibration was achieved by applying a shift to the zero-point 
from the 5577 A skyline. Cosmic rays were identified in the 
2-D spectra and edited out of the 1-D spectra. Finally, the 
red and blue spectra were joined together, averaging over 50 
A where the snr of the spectra were approximately equal. 

The final FWHM-aperture spectra are shown in Ap- 
pendix A. A gaussian fit was made to the lines to determine 
the line centres and FWHMs. Where there was significant 
absorption present in the line profile, this was removed be- 
fore the fit. The redshift of each quasar was determined from 
narrow line positions if possible, because these are more ac- 
curately measured and believed to be at the systemic red- 
shift. Where narrow line centres could not be measured ac- 
curately, a weighted mean of broad line centres was used to 
determine the quasar redshift. For some objects blue- ward 
absorption of Lya occurred, which shifts the measured line 
centre to higher wavelengths. In these cases the redshift mea- 
sured from Lya was not used in the redshift estimation pro- 
cedure. Line fluxes were measured with IRAF, by fitting a 
linear continuum and integrating the flux above this contin- 
uum. The observed equivalent widths were also measured in 
this way. The main source of error here is in accurately de- 
termining the underlying continuum. Appendix A contains 
the emission line data measured from the quasar spectra. 
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Name 


Telescope + 


Date 


Exposure 


olit wiatn 


ollt rTY 


1\T ^4- ^r-i 

IN otcs 


Previous 




Detector 




time (s) 


(arcsec) 


K ) 




reference 


5C6.5 


WHT+TSTS 

VV 11 _1_ | 1 l-J 1 l-J 


94Jan09 


600 


2.0 


75 


variable cloud cover 




5C6.8 




C/rtU dllVvO 


200 


2.0 


255 


vnriaKlo flrtnri rriupr 

V CLl 1CL L/1L LIU U LI LUVC1 






55 


95 Tan30 


1800 


2.0 


60 






5C6.34 


55 




500 


2.0 


75 


varianlo r'lonri pn\?f>T' 

V cLl IdUlL LIU LI LI LUVL1 




ICS 39 


„ 


97Tan09 


1200 


2.5 


141 






5C!fi q"i 


WHT4-FOS-2 

VY 11 1 |^1 v/U i 


91 Tanl 5 

£7 1 J till _1_ ij 


475 


4.0 


66 






5C!fi i fin 


WHT+TSTS 

VV 11 _L |^ J- 1 J J- i-J 


C/rtU cl 1 1 \J iy 


600 


2.0 


30 


V cLl IdUlL L1LJ L1L1 LLJ V LI 




^rjK 237 

•J \_< LI . iJiJ I 


55 


C/rtU CI 1 1 L/ O 


300 


2.0 


255 


varianlf 1 rlmin pn\?f>T' 

V CLl IdUlL L I LI LI LI L UVL1 


I 


5Cfi 251 


55 


94Jan08 


300 


3.0 


255 


variable cloud cover 




5C!fi 2fi4 


55 


C/rtU cl 1 1 \J iy 


1200 


2.0 


30 


varia nlo r"lriiiri rrn'or 

V CLl 1 cl [_/ 1 L L 1 LJ LI LI LUVL1 




5Cfi 282 


55 


95Jan31 


1300 


2.0 


81 






5C!fi 28fi 


55 




600 


2.0 


30 


V CLl IdL/lL L1LJL1L1 L UVL1 




5Cfi 287 


„ 




600 


2.0 


75 


ira r i a nl o plnnn pn\?f>T' 

V CLl IdL/lL LI LI LI LI LUVL1 




5Cfi 288 


55 


95Jan31 


1800 


2.7 


90 






5Cfi 291 


55 


97JanlO 


1800 


1.0 


24 






7C08084-2854 


„ 


97Anr07 


600 


2.5 


72 




2 


5C7.17 


55 


95 Tan29 

CsiJO f 11 1 _ ■ ' 


900 


3.1 


165 






5C7.70 


55 


97Apr07 


1200 


1.5 


120 






5C7.85 


55 


94JanlO 


730 


2.0 


280 


variable cloud cover 




5C7.87 


55 


95Jan30 


1700 


3.1 


60 






5C7.95 


LICK+KAST 


960ctl7 


600 


2.0 


143 






5C7.118 


WHT+ISIS 


97Apr07 


300 


2.5 


137 






5C7.194 


55 


97Apr07 


300 


2.5 


172 






5C7.195 


55 


95Feb01 


1200 


2.9 


70 






5C7.230 




97Apr07 


300 


2.5 


70 






7C0825+2930 




97Apr06 


900 


2.5 


170 







Table 2. Log of spectroscopic observations of the quasars. The ISIS spectra cover the range from about 3000 A to 8500 A with 
resolution approximately 10 A . The FOS-2 spectrograph at the WHT was used for 5C6.95. The spectrum of 5C7.95 was kindly obtained 
for us with the Kast spectrograph on the 3-m Shane Telescope at Lick Observatory by Andrew Bunker and Isobel Hook. Redshifts for 
5C6.237 and 7C0808+2854 had previously been determined in the literature . References: 1-Hook et al. (1996), 2-Barthel, Tytler & 
Thomson (1990). 



2.3 Optical photometry 

R and bj magnitudes have been measured for all the quasars 
and are listed in Table 3. If possible, these values were taken 
directly from the spectra. The FWZI aperture spectra were 
convolved with R and bj filter transmission profiles. The flux 
within the filter bandwidth was then measured. A similar 
procedure was followed with the spectrophotometry stan- 
dard, HD19445. This star has accurate photometry in the 
bj, R, I system (Gullixson et al. 1995). The measured quasar 
fluxes in each filter were then simply converted to a magni- 
tude. 

For the quasars observed in non-photometric condi- 
tions, this approach was not possible because a large fraction 
of the flux may be missing in the spectra. The variable cloud 
cover on these nights caused up to 2 magnitudes of extinc- 
tion in the spectra. However, the majority of these quasars 
are detected on the APM O and E plates and have measured 
magnitudes. In one case, the quasar was only detected on the 
O plate, so its R magnitude was estimated from its optical 
spectral index and bj magnitude. Two of the faint quasars 
were imaged in i?-band with the AUX port on the WHT, 
prior to spectroscopy. These were 5C6.33 and 5C6.39 which 
were observed on 1995 Jan 30 and 31, respectively. The ex- 
posure time for each image was 300s. A 300s J?-band image 
of 5C6.264 was obtained at the 107" telescope at the Mc- 



Donald Observatory, operated by the University of Texas 
at Austin, on 1995 Nov 20. These data were reduced us- 
ing standard procedures, as described in Lacy & Rawlings 
(1996). 

2.4 Galactic extinction corrections 

The 7C-I and 7C-II regions are at galactic latitudes, b ~ 
—30° and +30° , respectively. The amount of extinction and 
its variability across the sample was investigated. This is 
particularly important when one considers the intrinsic red- 
dening of quasars. IRA S 100/im maps of the regions were 
obtained from Skyview (http://skyview.gsfc.nasa.gov/) and 
the positions of the 7C radio sources over-plotted. The 
100/im surface brightness (corrected for zodiacal emission) 
in any direction is believed to be a very good indicator of 
the dust column density and therefore the visual extinction. 
The two variables are related by Av — 0.06 Zioo/im, where 
Aoo/im is in units of MJy sr _1 (Rowan- Robinson et al. 1991). 

The 100/xm maps show significant structure across each 
region with several clumps of galactic cirrus with peak fluxes 
more than twice the mean value. The mean extinction in 
each region is (A v ) = 0.24 (7G-I) and (Av) = 0.19 (7C- 
II). None of the quasars lie in the direction of the bright 
patches of cirrus. 5C6.95 is the quasar with the highest visual 
extinction value, Ay = 0.32. The difference between the B- 
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Table 3. Basic data and photometry for the quasar sample. Radio luminosities and absolute magnitudes were calculated assuming 
Qm = 1, = 0. K magnitudes are all measured using 5 arcsec diameter apertures except for 5C6.286, 5C6.288 & 7C0825+2930 (3, 
4 and 4 arcsec respectively, because of nearby objects), bj mags were calculated from B mags using B bj — 0.14 (e.g. Serjeant et 
al. 1998). Oopt, bj and R values corrected for galactic reddening: 7C-I, A B = 0.32, A R = 0.20; 7C-II, A B = 0.25, A R = 0.16. "The 
spectrum of 5C6.264 has a very poor snr because of several magnitudes of extinction by cloud and that of 5C6.288 only shows weak 
continuum because it is very faint. Therefore their optical spectral indices have been calculated assuming a power-law from their R and 
K mag nitudes. t5C7.17 & 5C7.118 are classified as broad-line radio galaxies (BLRGs), not quasars, because they have Mg > —23 (see 
Section 3). ^ These 5 quasars are core-jets sources (CJS; see Section 2). W5C7.230 is the only 7C quasar which does not have sufficient 
extended flux at 151 MHz to be included as a member of the complete sample. 
Key to ft j, R photometry: 

1 - APM catalogue 

2 - measured from spectrum 

3 - WHT AUX port R band image 

4 - McDonald Observatory, University of Texas at Austin, 107" Telescope R band image 

5 - used measured R mag. and assume bj — R = 0.5 adopting ct pt = 0.0 

6 - not detected on the E plate, only on the O plate, assumed bj — R = 0.5 to calculate R magnitude. 



band extinction to this quasar and the mean in 5C6 is 0.1 
mag. This is much less than typical magnitude errors (such 
as those from the APM), so we are justified in applying the 
same correction for objects in the same field. 

The effect of galactic extinction on the spectra was also 
investigated. A flat spectrum (a op t = 0) was created and 
then artificially reddened by various amounts (using the red- 
dening curve of Savage & Mathis 1979). It was found that 
there is a significant reddening of the initial flat spectrum. 
A power-law was fitted to the reddened spectra in the re- 
gion 4000-6000 A. For the least reddened spectra (those in 
7C-II), the reddened flat spectrum was fitted with a power- 
law with spectral index, a opt = 0.22. For 7C-I, a op t = 0.27 
and the most heavily galactic reddened quasar in the sam- 
ple, 5C6.95, a opt = 0.35. Hence, the scatter in a opt due to 



variable galactic reddening is < ±0.1 and therefore not sig- 
nificant compared with the errors i n th e determination of 
Qopt- The spectra shown in Figure A3 have not been de- 
reddened. However, the values of bj, R and a opt in Table 3 
have been corrected for galactic extinction. 



3 OPTICAL CLASSIFICATION OF 
RADIO-LOUD AGN 

Historically, quasars were defined as luminous objects with 
unresolved 'stellar' IDs and strong, broad (FWHM > 2000 
km s _1 ) emission lines. In contrast, radio galaxies have a 
resolved optical appearance and narrow (FWHM w 1000 
km s _1 ) or absent emission lines. However, it is now ap- 
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parent that this simple classification scheme is inadequate. 
The increased sensitivity and resolution of optical and near- 
IR imaging reveals host galaxies around quasars, particu- 
larly at low z. Hence a weak quasar in a luminous galaxy 
will have a resolved optical counterpart. Weak, broad emis- 
sion lines have also been seen in some radio galaxies. These 
objects are called broad-line radio galaxies (BLRGs) and 
have absolute magnitudes, Mb > —23. Spectropolarimetry 
of some radio galaxies has revealed broad lines scattered 
into our line-of-sight (e.g. 3C234: Antonucci 1984, Cygnus 
A: Ogle et al. 1997, 3C324: Cimatti et al. 1996), as expected 
in orientation-based unified schemes. It is also clear that 
some quasars have been reddened by dust so that the quasar 
is significantly dimmed in the UV (e.g. 3C22: Rawlings et 
al. 1995, 3C41: Simpson, Rawlings & Lacy, in prep.). The 
identification of this type of reddened quasar suffers from a 
redshift bias, because at higher redshifts optical spectra sam- 
ple shorter wavelengths where there is greater extinction. 
Near-infrared data is required to identify reddened quasars 
at high-z, e.g. the broad H-a line observed in the near-IR 
spectrum of 3C22 (Rawlings et al. 1995). Therefore the sim- 
ple requirement that quasars have broad emission lines must 
be treated with caution, because even with light reddening 
this becomes dependent on redshift and whether NIR data 
is available. 

Therefore we define quasars as objects with an unre- 
solved nuclear source which itself has Mb < — 23 (S1m = 1, 
Q.A = 0) after eliminating light from host galaxies and other 
diffuse continuum (such as radio-aligned emission, e.g. Mc- 
Carthy et al. 1987). This value is the same as the limit- 
ing absolute magnitude for quasars as defined in the quasar 
catalogue of Veron-Cetty & Veron (1996). Mb < —23 was 
decided upon after plotting Mb vs L151 for all 3CRR and 
7C Redshift Survey sources (see Section 5.1 for details of 
the Mb calculation). There is a fairly clear separation be- 
tween narrow-line only objects, which are less luminous than 
this, and objects with broad lines, which typically lie above 
Mb = —23. BLRGs fall below Mb = —23 and are conse- 
quently not classified as quasars. Further discussion of this 
classification scheme is deferred to Paper VIII. 

Two members of the 7C Redshift Survey are in the 
BLRG category: 5C7.17 and 5C7.118. They will not be con- 
sidered here as members of the complete quasar sample, al- 
though data is presented on them in Section 2. These objects 
are probably intrinsically weak AGN. They are not signif- 
icantly dust-reddened, optically-luminous quasars because 
their optical-UV continua are blue (see Appendix A). The 
only remaining problem with this definition is that heav- 
ily reddened quasars (Av S> 1) will be excluded because 
their optical luminosities are diminished to below our limit- 
ing absolute magnitude. Such objects certainly exist (Simp- 
son, Rawlings & Lacy, in prep.) and the consequences of this 
will be discussed in detail in Paper VIII. 

There are 2 sources in the 7C Redshift Survey which 
have optical luminosities Mb < —23 and which are spa- 
tially compact at about arcsec resolution, but for which we 
have no evidence of broad emission lines. We believe that 
both these objects are quasars and they are included as such 
in our sample. The spectrum of 5C6.264 was obtained in 
poor atmospheric conditions and we believe this is the rea- 
son that no broad lines are found in the spectrum. 5C6.288 
is at high redshift (z = 2.982) and shows only narrow lines 



in its optical spectrum. These objects are briefly described 
in Appendix A and will be discussed in more detail in Pa- 
per VIII. The optical spectrum of 5C7.195 also shows strong 
narrow lines with only a marginal broad component to the 
C III] emission line. This object is however known to be a 
lightly reddened (Av « 1) quasar (Willott et al. 1998a and 
Paper VIII). 

There are 8 objects traditionally classed as galaxies in 
the 3CRR catalogue which have integrated absolute mag- 
nitudes, Mb < —23. All of these are in the range —24 < 
Mb < —23 and subtraction of the host galaxy luminosity 
and/or extended emission provides a limit on their unre- 
solved optical component below our quasar definition value 
of Mb- The most luminous galaxy 3C226 has been shown 
by the HST imaging of Best, Longair & Rottgering (1998) 
to contain a maximum point-source contribution of 21% (at 
the 90% confidence level). Therefore the maximum absolute 
magnitude of its nucleus is —22.0, hence it is not by our 
definition a quasar. 

As mentioned above, near-infrared observations of 3C22 
and 3C41 (both listed by Spinrad et al. 1985 as radio 
galaxies) reveal reddened quasar nuclei. However, without 
de-reddening, these nuclei are fainter than our limit. Uni- 
fied schemes predict that all powerful radio galaxies con- 
tain quasar nuclei which are heavily reddened. Therefore we 
must exclude all reddened quasars with observed Mb > — 23 
nuclei from our sample. 3C109 has often been quoted as 
a BLRG because at z = 0.3056, its host galaxy is visi- 
ble in optical images. However, its nucleus is brighter than 
our limit of Mb = —23, so we classify this source as a 
quasar. Goodrich & Cohen (1992) estimate that there are 
~ 3 magnitudes of extinction towards the broad-line region 
of 3C109, so it is actually an intrinsically luminous quasar 
(Mb ~ —26). It should also be noted here that some RLQs 
show significant time variability in their nuclear optical lu- 
minosities and this may lead to some uncertainty in object 
classification. The full list of 3CRR quasars and BLRGs is 
in Appendix B. 



4 RLQ RADIO LUMINOSITY FUNCTION 

We can use complete radio-selected quasar samples to derive 
the radio luminosity function (RLF) of radio-loud quasars. 
The dataset here is a combination of the 40 quasars in the 
3CRR sample (Appendix B) and the 23 7C quasars with 
orientation-independent flux above the flux limit (as dis- 
cussed in Section 2). The 3 flat-spectrum quasars (3C345, 
3C454.3 & 5C7.230) which are promoted into the samples 
by their prominent cores have not been included. Whilst this 
is only a small sample (N = 63), the simplicity of the se- 
lection criteria compensate. The Si7g > 10.9 Jy flux limit 
of the 3CRR sample was recalculated at the same selection 
frequency as the 7C survey (assuming low-frequency spec- 
tral index, 0151 = 0.8), giving S151 > 12.43 Jy. The radio 
luminosities of the 3CRR and 7C sources were calculated 
at rest-frame 151 MHz using the spectral fitting technique 
described in Blundell et al. (1998a). The radio luminosity- 
redshift plane for both samples is shown in Figure nj. 
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Figure 1. The radio luminosity— redshift plane for the 3CRR and 7C quasar samples. 3CRR quasars are plotted as circles; filled are 
SSQs, open is the CJS. The triangles are 7C quasars; filled are SSQs, open are CJSs. Note that the combination of 3CRR and 7C samples 
removes the strong luminosity-redshift correlation present in a single flux-limited sample. The bins in luminosity and redshift used in 
Section 4.f are marked with dashed lines. Next to each bin is the mean redshift of the quasars in the bin. This plot is for cosmology I 

(n M = 1, n A = o). 



4.1 Binned 1/V a method 

The simplest way to estimate the luminosity function from a 
sample with known fluxes and redshifts is to use the binned 
1/Va method as described in Boyle et al. (1987). VI is the 
co-moving volume available for a source j in a bin with size 
Az (see Avni & Bahcall, 1980, for more details). p(L,z) is 
the number of sources per unit co-moving volume per unit 
log 10 Li5i. 

The RLF is given by 

JV 

i=l a 

From Marshall (1985), the error bars on each bin are given 
by 

/ N \ 05 

a =[j2( v «y 2 j (Aiog^isir 1 . 

To reduce the effects of small number statistics, we chose to 
have two redshift bins, equal in A log 10 (l + z). The redshift 
bins are 0.2 < z < 1.2 and 1.2 < z < 3.0. Each bin in 
luminosity is of size A log 10 L151 = 0.5. These bins are shown 
in Fig. @. 

The results of the binned analysis are shown in Fig. U 



(for cosmology I). For each redshift bin, the RLF appears to 
be consistent with a single power-law. The slope of the lu- 
minosity function, cti , appears to steepen considerably from 
Qi = 0.8 in the low-z bin to ai = 1.8 in the high- z bin. Mar- 
shall (1985) pointed out that with few redshift bins, there 
will be significant evolution across each bin. Since our cov- 
erage of the L151 — z plane is non-uniform, this leads to the 
mean redshift of objects in high-luminosity bins usually be- 
ing greater than that in the corresponding low- luminosity 
bins. Note there are exceptions to this, e.g. the two high- 
est luminosity, high- z bins contain only 3CRR quasars and 
consequently have lower mean redshifts than the lower lumi- 
nosity bins containing only 7C quasars. The peak in mean 
redshift in the bin with the highest redshift 7C quasars is 
the cause of the apparent bump in the binned high-z RLF. 
The mean redshifts of the sources in each bin are shown in 
Fig. |l| Due to the strong evolution from z ~ to z ~ 2 of 
RLQs, the main effect of this is an artificial flattening of the 
power-law slopes in the binned method. Hence the slopes 
quoted above are probably lower limits to the real values. 

For cosmology II, the binned luminosity function does 
not appear to be such a straight power-law as for cosmology 
I. The slope steepens slightly at high luminosities. This is 
most noticeable for the high-z bins (but may be partly an 
artefact of the binning method as discussed above). For cos- 
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mology III, the binned data has similar, straight slopes as 
for cosmology I. 

4.2 Parametric model fitting 

The binned analysis used in Section 4.1 is not good statisti- 
cal use of complete sample data because information is lost 
when the objects are binned. This is especially important 
when the most information possible is to be extracted from a 
small dataset. An analytic model of the luminosity function, 
p(L, z), is required. The model should have as few parame- 
ters as possible because of the small sample size being used 
to constrain it. The first simple model tested (model A) is 
a single power-law in luminosity with the redshift evolution 
specified by a gaussian, i.e. 




where p a is a normalising term, L = 10 27 W Hz _1 sr _1 , a\ 
is the power-law exponent, z a is the redshift of the gaussian 
peak and zi is the characteristic width of the gaussian. Thus 
there are only 4 free parameters in model A. The use of 
a gaussian redshift distribution is based on the observed 
redshift distribution of optically-selected quasars (Shaver et 
al. 1998) and the global star-formation rate (Madau et al. 
1996). 

To find the best-fit values of these parameters, the max- 
imum likelihood method of Marshall et al. (1983) was used. 
S is defined as —2 In £, where C is the likelihood function and 
terms independent of the model parameters can be dropped. 
In this method, the aim is to minimise the value of S, which 
is given by: 

S= -2^1n[p(L i ,« j )] + 2 JJ p(L,z)tl(L,z)^-dzdL 

where (dV/dz)dz is the differential co-moving volume ele- 
ment, fl(L, z) is the sky area available from the samples for 
these values of L and z and p(L, z) is the model distribution 
being tested. In the first term of this equation, the sum is 
over all the N sources in the combined sample. The second 
term is simply the integrand of the model being tested and 
should give w 2N for good fits. The integration range for the 
second term is < z < 5, 24 < log 10 L151 < 30. Integrating 
the RLF up to z = 10, instead of z = 5, in the maximum 
likelihood fitting made no difference to any of the models, 
because the co-moving volume available becomes very small 
at high redshift. Initial starting values of the free parame- 
ters were estimated from the binned RLF. The parameters 
were optimised via the downhill simplex method (Press et 
al. 1992) . Errors on the best-fit parameters were determined 
assuming a \ 2 distribution of AS (= S — 5 m in) (Lampton, 
Margon & Bowyer 1976). The error on each parameter was 
calculated by setting AS = 1 and allowing the other param- 
eters to float. The errors quoted here are 68% confidence 
levels, as discussed in Boyle et al. (1988). Table 4 contains 
the best-fit parameters of all the models considered. 

To assess the goodness-of-fit of a particular model the 
2-dimensional Kolmogorov-Smirnov (KS) test was applied. 
This multi-dimensional version of the KS test was first pro- 
posed by Peacock (1983) and the version used is that of 
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Figure 2. The radio luminosity function of RLQs derived here 
(for cosmology I). The points with la error bars are the binned 
RLF of Section 4.1. The solid lines are the maximum-likelihood 
best-fit model C and the dotted lines model H, plotted for several 
values of z. 



Press et al. (1992). It gives the probability, Pks, that the 
model distribution is a true representation of the dataset. 
Note that the KS test can only be used to reject models, 
not prove that they fit well. Indeed, for values of Pks ^ 0.2 
the probabilities may not be accurate and should not be 
used to distinguish between models. Therefore, we also eval- 
uate the likelihood directly from the minimisation routine 
and use these to compare well-fitting models. For model A, 
Pks = 0.03 for cosmology I and is even lower for cosmologies 
II & III, which indicates a poor fit. Another way to check 
the goodness-of-fit of a RLF model is to use a Monte-Carlo 
simulation to create an artificial L — z plane. These artificial 
datasets can then be compared to the actual L — z plane for 
the 3CRR and 7C quasar samples. For model A, the simu- 
lated data contained too many high radio luminosity sources 
(i.e. in 3CRR) at both low (z < 0.3) and high (2 < z < 3) 
redshift. 

An additional constraint on our models down to £151 « 
0.1 Jy comes from the 151 MHz source counts from the 7C 
and 6C surveys (McGilchrist et al. 1990 and Hales, Bald- 
win & Warner 1988, respectively). Since these counts include 
both quasars and radio galaxies, the source counts predicted 
from our quasar RLF models must be significantly lower. Al- 
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Figure 3. The normalised differential source counts from the 6C 
and 7C surveys are plotted here as solid lines and labelled 6C SC 
and 7C SC. Note that these counts include all radio sources, both 
quasars and radio galaxies. The solid circles with error bars show 
the binned quasar source counts from the 7C and 3CRR quasar 
samples described here. The source counts are normalised to the 
differential source counts for a uniform distribution in a Euclidean 
universe, such that dN a = 2400 (•S"^ 5 — Smii ) , where S" m i n and 
S max are the lower and upper flux limits of the bin. The other lines 
are quasar source count predictions from RLF models: A (short 
dash - long dash), B with log 1 g(Lbreak) = 26.6 (short dash), B 
with log 10 (L break ) = 26.8 (dot dash), B with log 10 (L broak ) = 
27.0 (long dash) and C with log 10 (L brcak ) = 26.8 (dot). 



though the fraction of quasars in complete samples at fluxes 
as low as 0.1 Jy is unknown, one would not expect it to be 
vastly different from at 0.5 Jy. Models are considered con- 
sistent with the source counts if the percentage of quasars 
predicted at 0.1 Jy is in the range 20%-40% (the quasar 
percentage in the 7C Redshift Survey (S151 > 0.5 Jy) is 
~ 30%). The single power-law model (A) predicts that the 
RLQ source counts are greater than the total 7C source 
counts at 0.1 Jy. Figure ^| plots the 6C and 7C total source 
counts and the quasar source counts predicted by several 
RLF models. Taken together with the poor goodness-of-fit 
test from the L — z data, the source counts constraint forces 
us to conclude that model A should be rejected. 

The steep-spectrum RLF of Dunlop & Peacock (1990) 
(which includes radio galaxies and quasars) could be fitted 

i L2.7 ~ 



25.5. This is equivalent to log 10 L151 ta 26.5 (assuming 
arad = 0.8). The over-prediction of the RLQ source counts 
by model A could be caused by the lack of a break in 
the luminosity function. This is because samples selected 
with 5*151 > 0.1 Jy, would contain low luminosity sources 



(log 10 L151 ta 26) out to redshifts, 



1.5. Therefore, the 



excess counts are due to the increased number density of 
low radio luminosity sources in models with an unbroken 
power- law. 

Therefore, our next model (B) is a broken power-law 
in luminosity with a break at Lbreak. Unfortunately, we are 
unable to derive the value of Lbrcak or the power-law index 
below the break with our complete sample data, because we 
have so few quasars at these low radio luminosities. We as- 
sume that the RLF is flat (ai = 0) below the break. The 
value of Lbrcak can be estimated using the 6C/7C source 
counts, because the counts at 0.1 Jy are very sensitive to 
it. For cosmology I we can constrain Lbrcak for this model 
to 26.7 < log 10 (Lbrcak) < 26.9 (see Fig. |). The quasar per- 
centages for log 10 (L b rcak) = 26.6 and log 10 (L br cak) = 27.0 
are 45% and 15%, respectively. Log 10 (Lbrcak) = 26.8 was 
assumed since it gives a quasar fraction at 0.1 Jy sim- 
ilar to that in the 7C sample considered here. Similarly 
log 10 (Lbreak) = 27.2 and 27.3 for cosmologies II and III, 
respectively. Note that the value of log 10 (Lbreak) is not a 
free parameter in our fits - it is constrained by the source 
counts, not the L151 — 2 data. Maximum likelihood fits to 
the 4 free parameters of model B gave values of Pks in the 
range 0.4-0.7 for all cosmologies which indicates that the 
model fits our data fairly well. 

The functional form of the redshift evolution of the RLF 
in these models is symmetric in z about z a . This form was 
not based on theoretical predictions, but merely because it 
appears to be a good fit to the increase in density from z = 
to 2 ~ 2 (although clearly p must decline at some high value 
of z). In the fit using model B above, z a = 2.2 (for cosmology 
I). There are only 6 quasars in the 3CRR/7C sample with a 
redshift greater than this. Therefore, we cannot claim that 
there is a symmetric decline in the space density at 2 > z , 
just because we get a high probability from our goodness- 
of-fit test. To test whether a high redshift decline in p is 
necessary, we repeated the model-fitting procedure using a 
one-tailed gaussian in 2, such that p(z) does not decline at 
all in the region z < 2 < 5 (and p — beyond 2 = 5, as 
before). This new model (C) fits the data equally well with 
similar KS probabilities and a slightly greater likelihood. 
The peak redshift for C ranges from 1.7 for cosmology I 
to 2.0 for cosmology III. Hence we have no evidence for a 
decline in p at 2 ~ 2 and we adopt model C as a slightly 
better fit to the data than B. 

Many derivations of the OLF and XLF have used the 
functional form (1 + z) k to describe the redshift evolution of 
the quasar population (e.g. Boyle et al. 1994; Goldschmidt 
& Miller 1998). We have replaced the gaussian redshift term 
from model A with this new term to make model D. Note 
that this reduces the number of free parameters in our fit 
from 4 to 3. Model D provides a very poor fit. The KS proba- 
bilities are lower than 10 _s for all cosmologies. Monte-Carlo 
simulations of this model create many sources at low redshift 
and high luminosities, completely unlike the actual L151 — 2 
distribution. As before, we now include a low-luminosity 
break in the RLF (model E, analogous to model B). This 
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Errors for model C (cosmology I): log 10 (p ); -0.12, +0.13, «i ; -0.09, +0.09, z c ; -0.19, +0.22, zi ; -0.06, +0.07. 

Table 4. Best-fit parameters for model RLQ RLFs, as described in Section 4.2. Cos is the cosmological model assumed, 7V par is the 
number of free parameters in the fit, C is the likelihood from the fit, relative to model C for each cosmology and Q% is the quasar 
percentage predicted at 0.1 Jy from the source counts. A brief description of each model is given below (see Section 4.2 for more details). 
A: single power-law RLF, 2-tailcd gaussian z dependence. 

B: single power-law RLF with low luminosity break, 2-tailed gaussian z dependence. 
C: single power-law RLF with low luminosity break, 1-tailed gaussian z dependence. 
D: single power-law RLF, (1 + z) k z-dcpcndcncc. 

E: single power-law RLF with low luminosity break, (1 + z) k z dependence. 
F: single power-law RLF with low luminosity break, (1 + z) k z dependence with z max = z . 
G: z dependent power-law RLF with low luminosity break, 2-tailed gaussian z dependence. 
H: z dependent power-law RLF with low luminosity break, 1-tailcd gaussian z dependence. 



model also fits poorly with Pks < 0.05 for all cosmologies. 
The relative likelihoods between model E and model C are 
< 10~ 6 . The low-luminosity break in model E is required to 
be 0.4 dex greater than models B & C and even then the per- 
centage of quasars in the source counts at 0.1 Jy are > 50%. 
Hence model E must be rejected. The reason for its poor fit 
is that the number density of quasars continues increasing 
up to high redshifts. This is illustrated by the Monte-Carlo 
simulations which show many sources at z > 3, unlike the 
real data. 

Boyle et al. (1994) find a good fit to the XLF with 
(1 + z) k evolution if one allows a maximum redshift for 
this evolution, z a , such that p is constant at z > z a . This 
should solve the problem of model E - too many high red- 
shift quasars. Model F implements this form for the evolu- 
tion and consequently there are now 4 free parameters in 
the fit. Model F is now extremely similar to model C, in 
that in both p increases up to z and then remains con- 
stant. The fits for model F do indeed have similar parame- 
ter values to model C. The KS probabilities are in the range 
0.6 — 0.8. Model F agrees well with the source counts (with 
log 10 (Lbrcak) = 26.8, as for model C) and Monte-Carlo sim- 



ulations. The relative likelihoods of model F to model C 
are a factor of 2 to 5 lower, so wc retain model C as the 
best-fitting model. 

Our binned RLQ RLF has a power-law in L151 which 
apparently steepens as the redshift increases. Therefore, the 
next logical step is to include this possibility into our RLF 
(model G). The power-law part of the model RLF is now 

for L151 > Lbroak and a flat power-law (ai = 0) for 
£151 < ibrcak- This adds another free parameter to our 
model as Qi has been replaced by a\ + a^z. However, 5 
free parameters is still a small number considering this is a 
3-dimensional model: p (L,z). The redshift distribution is 
again specified by a two-tailed gaussian, as in model B. The 
2-D KS test was applied to the best-fits to this model and 
gave Pks = 0.4 — 0.7 for all the cosmological models. Monte- 
Carlo simulations using the model G RLF give L — z planes 
indistinguishable from the true data. In addition, this model 
also fits the source counts at low fluxes well, giving a quasar 
fraction at 0.1 Jy of w 30%. Similarly, we tested a variant 
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Figure 4. The plot at top-left shows the L151 — z plane for the 3CRR and 7C quasar samples. The other 3 plots show Monte-Carlo 
simulations of the L151 — z plane for samples with the same flux limits and sky areas as the 3CRR and 7C samples. These simulations 
were generated using the best-fit RLF models B, C and H. All plots are for cosmology I. 



of model G with a one-tailed gaussian redshift distribution 
(i.e. no decline in redshift for z > z a ). This model (H) again 
fits well with KS probabilities in the range 0.4 — 0.9 and 
good Monte-Carlo simulations and source counts fits. The 
likelihoods for model H are slightly higher than for model 
G. 

Comparing the relative likelihoods for models C and H, 
we find that for cosmology I model H has a 10 times greater 
likelihood than C, whereas for cosmologies II and III the 
likelihoods are comparable. Since H is a 5 parameter model, 
whereas C has only 4 parameters, we should only accept H 
as a best-fitting model if it can be shown that it is a sig- 
nificantly better fit than C. To do this we need the ratio of 
the posterior probabilities for each model confronted with 
the data. The required ratio can be fairly approximated by 
the ratio of the likelihoods multiplied by an 'Ockham factor' 
(< 1) which penalises the extra latitude afforded to model 
H by an extra free parameter (e.g. Sivia 1996). For cosmolo- 
gies II and III there is clearly no motivation to prefer model 
H (since their likelihood ratios are about 1). For cosmology 
I the ratio of the posterior probabilities is at most 10, so 
although this hints that model H is to be preferred, uncer- 
tainty in the value of the Ockham factor ensures this remains 
a hint rather than a firm conclusion. The simulated L151 — z 
data and source counts plots for both models are similar and 
fit the data well. Note that there are no 3CRR RLQs with 
z > 2.02, despite the fact that the highest redshift RLQ, 



3C9, is luminous enough to be seen in the 3CRR sample at 
z > 3 and 3C191 at z « 2.7. In contrast, there are nine 
7C RLQs in the range 2 < z < 3. This difference is repro- 
duced slightly better by model H than model C, due to the 
steepening of the RLF at high-z. An additional check on the 
validity of introducing an extra parameter comes from the 
size of its error. For model H, the extra parameter has the 
value 0L2 = 0.46lg'29- These errors are approximately equiv- 
alent to la errors, so for the additional term to become zero 
(and hence revert to a 4 parameter model) would only be a 
1.6cr result. Therefore, we conclude that we do not have suf- 
ficient evidence for the 5 parameter models to be considered 
a better fit than the 4 parameter models and hence model 
C remains our preferred model until further data suggest 
otherwise. 

4.3 Results 

The RLF models C and H are plotted in Figure ^ for z = 0.5, 
1, 2 and 3 (for cosmology I), along with the binned RLF of 
Section 4.1. Although at first glance, the binned and model 
RLFs appear quite different, these apparent differences are 
easily explained by the fact that the mean redshift varies 
from bin to bin, as discussed in Section 4.1. The fairly flat 
slope of the low- z binned RLF is due to the mean redshift 
increasing from 0.5 to TO as luminosity increases. This can 
clearly be seen by comparing the binned points with the 
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model lines for z = 0.5 and z = 1.0. At high redshift (1.2 < 
z < 3.0), both model C and H lines slightly underestimate 
the number density in several of the bins. This may simply 
be a consequence of small number statistics, because the 
largest deviation is only ~ 2a. The small sample size used 
here means that the RLF is only well-constrained in regions 
of the Lisi — z plane populated by 3CRR and 7C sources 
and extrapolations to other ranges of luminosity and redshift 
should be regarded with caution. 

It is however, useful to consider extrapolation of the 
RLF to lower radio luminosities than those found in this 
sample. It has long been known that the radio luminosity 
distribution of quasars is bi-modal, with the separation into 
radio-quiet and radio-loud quasars occurring at log 10 (Lisi) 
~ 25.5 (c.f. Kellermann et al. 1989). The 6C and 7C source 
counts show that a break in the power-law is required at 
log 10 (Li5i) » 27, although we cannot constrain the slope of 
the RLQ luminosity function on the low luminosity side of 
this break or the exact luminosity of the break. The paucity 
of quasars (only 2) in the combined 3CRR/7C sample with 
log 10 (Lisi) < 26.5 is consistent with a break in the RLF 
near this luminosity. It should be noted however that 7 of 
the BLRGs in 3CRR & 7C have log 10 (Li 5 i) < 26.5. In fact, 
due to the radio-optical correlation to be discussed in Sec- 
tion 5, BLRGs often have low radio luminosities and may 
contribute significantly to the quasar RLF at these radio lu- 
minosities. However, even including low optical luminosity 
quasars (BLRGs), there must still be a break in the RLQ 
RLF to be consistent with the 6C/7C source counts. 

The 6C complete sample of Eales ( 1985a, b) was selected 
at 151 MHz with flux limits 2.00 < Swi < 3.93 Jy and covers 
0.102 sr (using the revised selection criteria from Blundell et 
al. 1998a). Hence it lies in between the 3CRR and 7C sam- 
ples in the L151 — z plane and can be used to check our RLF 
models. It contains only seven quasars (Eales 1985b; unpub- 
lished data) , none of which are promoted into the sample by 
their core fluxes (Blundell et al. 1998a; unpublished data) . 
Integrating our model C RLF (for cosmology I) over the 6C 
flux limits and sky area, we find that there should be, on 
average, 16 quasars in the sample. Assuming Poisson statis- 
tics, p(< 7) = 0.01 given mean /1 = 16. The result is similar 
for model H with 15 quasars predicted in the 6C sample. It 
is therefore just plausible that the low quasar fraction in this 
sample (as commented on by Eales) is simply a statistical 
fluke. 

The radio-loud quasar RLF derived here has some re- 
markably similar features to the OLF and XLF. The in- 
crease in log 10 p(z) from 2 ~ 0.5 to z ~ 2 for model C is 
~ 1.5, compared to w 1.8 for medium luminosity quasars 
in the XLF (Boyle et al. 1994) and sa 1.5 for medium lu- 
minosity quasars in the OLF (Goldschmidt & Miller 1998). 
Boyle et al. (1994) fit the quasar XLF with a redshift dis- 
tribution which peaks at 2 = 1.6 and then remains constant 
at higher redshifts, in accordance with our best-fitting RLF 
models. Of course, a decline in the OLF beyond 2 m 3 is now 
established (e.g. Warren, Hewett & Osmer 1994; Schmidt, 
Schneider & Gunn 1995), although these optically-selected 
samples could be biased if there is significant dust at high 
redshift. Flat-spectrum radio quasars also show decreasing 
number densities beyond 2 ~ 3 (Shaver et al. 1998), a result 
which is insensitive to the effects of dust. The high lumi- 
nosity part of the OLF is known to steepen at high redshift 



(Goldschmidt & Miller 1998) and we find weak evidence 
that this may be occurring also in the RLF (e.g. model H). 
These similarities are important to pursue further because 
radio sources are known to be short-lived compared to the 
Hubble time (typically 10 7 yr; Scheuer 1995), so the evo- 
lution must really be telling us about the evolution of the 
birth-rate (see also Blundell et al. 1998a), and by analogy 
the same for the OLF. 



5 THE RADIO-OPTICAL CORRELATION 

5.1 Absolute B magnitude calculation 

iC-corrections in the B-band for high redshift objects can 
be very large and uncertain. The multi-colour photometry 
presented here avoids this potential source of error, because 
it allows the luminosity of an object in the rest-frame B- 
band to be calculated. For example, at 2 ~ 3 rest-frame 
B-band emission is observed at 1.8//m (close to the near-IR 
K-ba,nd). For sources at 2 > 0.5, emitted _B-band radiation 
is observed between the R and K bands. Therefore, for all 
7C quasars a power-law continuum was fitted between the R 
and K magnitudes and a rest-frame B-band luminosity cal- 
culated. This was then converted to an absolute magnitude, 
Mb- The errors on these values are typically ±0.2 mag; the 
main source of error being the B and R magnitudes from 
the APM catalogue. Paper VIII shows that this method pro- 
vides a good estimate of Mb, even in quasars reddened by 
A v « 1. 

For the 3CRR quasars, complete A"-band photometry 
has not been published, so this approach was not possible. 
However, the maximum redshift in 3CRR is ~ 2, so the K- 
correction terms are not so large. For the 3CRR quasars a 
power-law optical continuum was calculated using published 
V and R magnitudes or V and B magnitudes, where avail- 
able. For the few sources with only a published V magnitude, 
an optical continuum slope of a opt = 0.5 was assumed. This 
slope is equivalent to B — K ~ 3 and is typical of unred- 
dened quasars (Francis 1996) . The calculated absolute mag- 
nitudes are listed in Table Bl in Appendix B. Jf-magnitudes 
were used (where available) in the calculation of the absolute 
magnitudes of 3CRR radio galaxies referred to in Section 2. 

5.2 Results 

All the numerical results in this Section are presented for 
cosmology I. None of the conclusions differ significantly for 
either cosmology II or III. In Figure ^, we plot the cal- 
culated absolute magnitudes of the quasars and broad-line 
radio galaxies against redshift. The horizontal line is the 
dividing line in absolute magnitude between quasars and 
BLRGs from Section 3. The 5 CJSs are considered sepa- 
rately from the SSQs in this section because of possible ori- 
entation biases in their optical luminosities (e.g. Browne & 
Wright 1985). The BLRGs are also excluded because they 
are probably not a complete sample. There is a very weak 
correlation between Mb and 2 for the SSQs (in the sense 
that the quasars increase in luminosity at greater 2). Using 
Spearman's method of rank correlation, the significance of 
the correlation between Mb and 2 is 88%. A radio-optical 
correlation is implicit in Figure H in the redshift range 
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Figure 5. Absolute B magnitude against redshift for the 3CRR 
and 7C quasar samples. The plotting symbols are the same as 
those used in Figure [j]. BLRGs with z > 0.1 are also plotted here 
(asterisks - 3CRR, plus signs - 7C). This plot is for cosmology I. 
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Figure 6. Absolute B magnitude against low-frequency radio 
luminosity for the 3CRR and 7C quasar samples. The plotting 
symbols are the same as those used in Figure |i| Plotted also is 
the best-fit straight line to the SSQs only, as described in the 
text (solid line). The dashed line is the radio-optical correlation 
for SSQs of Serjeant et al. (1998) transformed to cosmology I and 
Ho = 50 km s _1 Mpc -1 . This plot is for cosmology I. 



(1 < z < 2) where there is the most overlap between the 
3CRR and 7C samples. It can clearly be seen that the 3CRR 
quasars (which have radio luminosities ~ 25 times greater 
than 7C quasars at the same z) are systematically brighter 
than those in 7C. In this redshift range, the mean absolute 
magnitude, (Mb), for 3CRR SSQs is —25.8 with standard 
error 0.2 mag. The 7C SSQs have (Ms) = -24.6 ±0.3 mag, 
so are significantly fainter by about one magnitude. Note 
that the BLRGs appear to be more tightly correlated than 
the quasars on this plot. This is most likely due to the fact 
that all but 2 BLRGs are from the 3CRR sample, so there 
is a tight Z/151 — z correlation. 

Figure ^| shows the correlation between radio and opti- 
cal luminosities for the combined 3CRR and 7C SSQ sam- 



ples. It can clearly be seen that there is a much stronger 
correlation than in Figure ^. The significance of this corre- 
lation is > 99.9%. This correlation does not appear to be 
caused by a correlation between Mb and z (and a L\ \ — 2 
correlation), because the statistical significance is very much 
greater. Also plotted on Figure |(S| is the best-fit straight line 
to the 3CRR/7C SSQs (calculated by minimising the sum 
of the squares of the residuals). This line has a slope of 
—2.5 x (0.4 ± 0.1), i.e. L opt oc £rad- The dispersion about 
our straight-line fit is 1.0 magnitude. This is in fairly good 
agreement with Serjeant et al. (1998) who found a slightly 
steeper slope of —2.5 x (0.6 ± 0.1), which is also plotted 
in Fig. ^ for comparison. One possible reason for our find- 
ing a flatter slope than Serjeant et al. follows from their 
lack of a limiting absolute magnitude in their definition of 
a quasar. Their sample contains quasars with Mb > —23 
(which we would classify as BLRGs). These objects mostly 
have low radio luminosities, therefore increasing the slope 
of the correlation. This is evident in Fig. |(| There may be 
other systematic differences between the RLQs in the Ser- 
jeant et al. MAQS sample and those in the 7C sample. For 
example, the MAQS sample is selected at 408 MHz lead- 
ing to a higher proportion of compact steep-spectrum radio 
sources. We will consider these points in more detail in a 
future paper. 

To investigate whether the radio-optical correlation 
could be due to a Mb — z correlation, the Spearman par- 
tial rank correlation coefficient was used. This assesses the 
statistical significance of correlations between two variables 
in the presence of a third (Macklin 1982). For the correla- 
tion between Mb and 2 (independent of £151), tm b z,l 151 = 
—0.13 with significance Dm b z,l 151 = —1.0. The significance 
is equivalent to the deviation from a unit variance normal 
distribution if there is no correlation present. Hence the cor- 
relation between Mb and 2 is marginally significant. Now we 
consider the correlation between Mb and L151 (independent 
of z). We find r MB L 151 ,z = -0.43 and D Mb l 151 ,z = -3.6. 
Hence we can be sure that the obvious radio-optical corre- 
lation of Figure ^| is not caused by the Mb — z correlation 
which plagues single flux-limited samples. The high radio 
luminosity but low optical luminosity (but brighter than 
Mb = —23) quasar region of Figure ^|is genuinely devoid of 
quasars. There are no selection effects which would cause a 
lack of quasars in this region. 

Whenever a correlation such as this is presented, it is 
essential that all possible sources of error, both random and 
systematic, are considered. This is particularly important 
when combining samples with different selection criteria. 
The 3CRR and 7C SSQ samples are however, selected with 
very similar criteria (apart from their different flux limits). 
The only difference is their slightly different selection fre- 
quencies; 3CRR is selected at 178 MHz and 7C at 151 MHz. 
Both are optically complete, because all the sources in the 
3CRR sample and the 7C Redshift Survey have been identi- 
fied. Hence we have eliminated one of the previously impor- 
tant possible sources of error, that of 'missing' quasars due 
to an optical magnitude limit. 

There may be systematic errors between the two sam- 
ples because the optical data available for them is inhomo- 
geneous. The right side of Figure ^| is composed mostly of 
3CRR objects and the left side mostly 7C, so could this ex- 
plain the observed correlation? The absolute magnitudes of 
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the 7C quasars were calculated using multi-colour photom- 
etry, so are likely to be good estimates of the true values. 
Errors in Mb due to reddening are not expected to be cru- 
cial here, because R and K magnitudes were predominantly 
used. Intrinsic reddening can be important when using B 
magnitudes for high redshift objects, because the observed 
radiation was emitted in the UV. For the 3CRR quasars, Mb 
was calculated using optical photometry. Hence it is possible 
here that some quasars may be significantly reddened and 
this is a possible source of error. In addition, the magnitudes 
listed in LRL have not been corrected for galactic extinction. 
However, both these effects would tend to cause our calcu- 
lated 3CRR Mb values to be fainter than their true values. 
This cannot explain the observed radio-optical correlation, 
because the 3CRR quasars are systematically brighter than 
those in 7C. As noted above, the mean absolute magnitude 
of 3CRR SSQs is about 1 magnitude greater than the mean 
for 7C SSQs in the same redshift range. Hence the 3CRR 
photometry must be systematically in error by this amount 
for the correlation to disappear. It seems unlikely that this 
could be the case. The numerical results presented in this 
section are virtually independent of the choice of cosmology. 

The data we have here is not detailed enough for us to 
explore the nature of the radio-optical correlation in quasars. 
For a discussion of possible physical causes of this effect, 
we refer the reader to Serjeant et al. (1998) and Paper III. 
Note that the CJSs tend to occupy positions on the optically 
bright side of the correlation. This could be explained by 
orientation-dependent optical emission, however the small 
number (6) of CJSs in the combined sample prevents us 
from drawing any firm conclusions. 



6 COMPLETENESS OF OTHER SAMPLES OF 
SSQS 

Using the quasar radio luminosity function derived in the 
previous section (model C), artificial RLQ samples were gen- 
erated with a Monte-Carlo simulation for a range of S151 flux 
limits from 0.1 to 20 Jy. Each sample contains approximately 
5 x 10 4 sources, each defined by L151 and z. The straight- 
line fit to the radio-optical correlation for SSQs derived in 
Section 5 defines a characteristic value of Mb for each value 
of L151 . The scatter about this straight-line was assumed to 
be independent of radio luminosity and a gaussian scatter of 
1.0 magnitudes introduced into the Mb calculation. For each 
source, Mb and z were then used to calculate its apparent B 
magnitude (assuming a opt = 0.5). Quasars with Mb > —23 
(BLRGs) were excluded from the simulated samples. The 
fraction of quasars brighter than the limiting B magnitude 
Slim was then calculated for each sample as a function of 

Slim- 

Figure |^ shows contours of the percentage of SSQs ex- 
pected to be detected in surveys as a function of their optical 
and radio flux limits (for Qm = 1, Qa = 0). The immediate 
conclusion to be drawn here is that as surveys go to fainter 
radio fluxes, they must also go deeper optically, otherwise 
a significant fraction of quasars are missed. For Qm — 0, 
Qa = 0, the contours have steeper slopes across the diagram, 
such that they diverge from the Qm = 1, Qa = contours 
at faint radio fluxes to be 0.5 magnitudes in B fainter at 
S151 = 0.1 Jy. A similar behaviour is seen for Qm = 0.1, 



Qa = 0.9. It must be remembered that these contours have 
been calculated assuming models for the radio-optical cor- 
relation and the quasar RLF from the combined 3CRR/7C 
sample. Therefore in regions not probed by these samples, it 
is an extrapolation of the models. The area with the largest 
uncertainty is therefore S151 < 0.5 Jy. At S151 =0.1 Jy, 
typical errors on the contours may be as large as ±1.0 in 
limiting B magnitude. 

According to Figure |, 99% of the 3CRR SSQ sample 
has B < 20.5 (i.e. brighter than the POSS limit). This is in 
good agreement with the fact that all of them are detected 
on the POSS plates. Our calculation predicts that 13 out 
of the 18 7C steep-spectrum quasars have B < 20.5, com- 
pared with the actual number which is 11. Note that the 
completeness percentages calculated here are strictly upper 
limits, because they do not take fully into account the faint- 
ness of reddened quasars. The completeness contours de- 
rived here only apply to low-frequency selected samples, e.g. 
151 MHz. This is because high-frequency samples contain a 
large fraction of FSQs which are typically brighter than the 
mean SSQ radio-optical correlation (Wills & Lynds, 1978; 
Browne & Wright, 1985). Therefore FSQs are more likely to 
be brighter than the optical magnitude limits of samples. 

We can use Figure ^ to estimate the completeness of 
various low-frequency radio selected quasar samples. This 
allows us to investigate possible bias in the sample proper- 
ties due to selection effects. The large 7C quasar sample of 
Riley et al. (1998) has a faint radio flux limit of S151 > 0.1 
Jy and a bright optical limit of 7? < 20. This sample is prob- 
ably only approximately 40% complete. There are other new 
radio quasar samples which use POSS-I (R m 20) identifica- 
tions. The B3-VLA quasar sample (Vigotti et al. 1997) with 
a flux limit S408 > 0.1 Jy is estimated to be 55% complete 
from our plot (c.f. their estimate based on optical magnitude 
histograms of 70%). For samples with S151 ^ 1 Jy, one must 
use optical data significantly deeper than the POSS-I limit 
in order to avoid severe incompleteness in these samples. 

The MAQS sample of Serjeant et al. (1998) is selected 
by S408 > 0.95 Jy and B < 22.5. We estimate that this sam- 
ple is 99% complete, due to its faint optical limit (and its 
relatively bright radio limit). However, the higher selection 
frequency of 408 MHz samples rest-frame > 1 GHz emission 
for z > 1.5 quasars. At these rest-frame frequencies cores 
and jets often dominate over extended emission, so the ex- 
tended radio luminosity is lower and hence the quasars will 
be lower down the radio-optical correlation. It is possible 
however that a beamed optical-UV component, as discussed 
in the previous paragraph, may compensate in these objects. 
Therefore we conclude that the MAQS sample is likely to 
have excluded very few quasars due to its optical magnitude 
limit. A comparison of our RLF models with the MAQS sam- 
ple will appear elsewhere (Serjeant et al. in preparation). 



7 CONCLUSIONS AND SUMMARY 

The sample of 7C quasars presented here is complete because 
every source in the 7C Redshift Survey has been identified. 
Therefore we do not suffer from colour biases excluding red 
and/or faint quasars. However, if orientation-based unified 
schemes (e.g. Antonucci 1993) are correct, then powerful ra- 
dio galaxies are reddened quasars. For intermediate values 
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Figure 7. Percentage of quasars detected in surveys as a function of their optical and radio flux limits. This plot is for cosmology I. 



of reddening (Av ~ 2), different observers may classify ob- 
jects as either radio galaxies or quasars. These issues will be 
investigated further in Paper VIII. The main conclusions to 
be drawn from this paper are: 



or there will be serious incompleteness in the samples. This 
would introduce many biases into the analysis of such sam- 
ples, because only the optically brightest quasars would be 
included. 



• The radio luminosity function of radio-loud quasars has 
been derived. The best-fitting model is a single power-law 
with slope tvi = 1.9 ± 0.1. We find that there must be a 
break in the RLQ RLF at log 10 (Li 5 i/ W Hz" 1 sr" 1 ) <27, 
in order for the models to be consistent with the 7C and 
6C source counts. The z-dependence of the RLF follows a 
one-tailed gaussian which peaks at z = 1.7±0.2 and remains 
constant at higher redshifts. Models with or without a red- 
shift cut-off beyond z « 2 fit almost equally well. We find 
marginal evidence that the RLF slope steepens as redshift 
increases. This form for the evolving RLF is very similar to 
the OLF determined from optically-selected samples (e.g. 
Goldschmidt & Miller 1998). 

• We have confirmed the existence of a significant cor- 
relation between the low-frequency radio and optical lumi- 
nosities of SSQs, as reported in Serjeant et al. (1998). This 
correlation is not due to correlations with redshift or sample 
selection effects. 

• Using the derived radio-loud quasar RLF and the ob- 
served radio-optical correlation, we have been able to es- 
timate the fraction of quasars identified in low-frequency 
radio-selected samples with various magnitude limits. It is 
found that samples with faint radio flux limits must also go 
much deeper in the optical than the POSS-I plates (R « 20), 
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APPENDIX A: OPTICAL SPECTRA 



In this appendix we present the optical spectra and emis- 
sion line data for the 7C Redshift Su rvey quasars and broad- 
lined radio galaxies. Figures Al and A2 are the optical spec- 



tra of 5C6.264 and 5C6.288, respectively. These are the two 
quasars with no evidence of broad emission line s which will 
be discussed further in Paper VIII. Figure A3 displays the 



optical spectra for the other 24 7C quasars and broad-line 
radio galaxies. Prominent emission lines are labelled. Table 
Al lists the properties of the emission lines measured from 
all these spectra. 
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name 


z 


line 


Arest 


-^obs 


FWHM 


flux 


snr 


w A 




notes 








(A) 


(A) 


(A) 


(W m- 2 ) 


<T 


(A) 






5C6.5 


1.038 


C III] 


1909 


3887 ± 15 


68 ± 10 


1.2c-17 


5 


71 ± 


20 


abs 


lO Rn9nzL1A 44-^39991 


-L 

ZJZ 


P TT1 


2326 


/l 7/1 q 4- i n 


jo in -Lz 


9 1 fi 
Z.Uc-lo 


2 1 


1 Q -t- 


r: 
o 






.002 


Mg II 


2799 


5704 ± 2 


139 ± 40 


1.3e-17 


4 


138 ± 


40 




5C6.8 


1.213 


C IV 


1549 


3434 ± 7 


34 ± 20 


1.6c-17 


2.8 


65 ± 


35 


in noise 


7C B020439.6+313752 


± 


He II 


1640 


3626 ± 3 


23 ± 6 


3.4c-18 


1.6 


16 ± 


8 






nn9 


r* TTTl 


1 QflQ 




iin -t- 
iiu m ou 


1 4e-17 


K 
O 


£1 -t- 










Mg II 


2799 


6196 ± 5 


101 ± 20 


l.Oc-17 


6 


71 ± 


16 


abs 






[O II] 


3727 


8247 ± 5 


24 ± 5 


l.Oc-18 


1.3 


23 ± 


8 




5C6.33 


1.496 


C IV 


1549 


3855 ± 10 


150 ± 40 


1.9e-18 


7 


220 it 


70 


abs 


7C* Rfl9nS^fi 14-mzL1 Q 


4- 


He II 


J-O^U 


AOQ« -t- 1 O 


iin -|_ in 
iiu m ou 


4 Oc-19 


9 °, 
z.o 


^ -t- 
jo m 


1 7 






.003 


C III] 


1909 


4752 ± 20 


170 ± 40 


9.5e-19 


4 


229 it 


70 


abs 






Mg II 


2799 


6998 ± 20 


209 ± 50 


8.8e-19 


4 


200 it 


60 


abs 


5C6.34 


2.118 


Ly a 1 N V 


1216 


3790 ± 6 


47 ± 10 


5.3c-17 


35 


198 ± 


40 


abs 


7C B020838. 5+312131 


± 


C IV 


1549 


4827 ± 3 


76 ± 10 


1.5e-17 


15 


71 it 


10 


abs 




.002 


He II 


1640 


5116 ± 4 


42 ± 8 


6.1e-18 


6 


38 it 


12 








C III] 


1909 


5936 ± 15 


146 ± 20 


8.1c-18 


6 


56 ± 


30 


cr in line 


5C6.39 


1.437 


Si IV / O IV 


1400 


3408 ± 15 


57 ± 20 


8.9e-19 


7 


193 ± 120 




7C B020853.8+321446 


± 


C IV 


1549 


3774 ± 3 


65 ± 12 


6.1c-18 


70 


415 ± 


80 






.003 


He II 


1640 


3994 ± 3 


17 ± 7 


3.5c-19 


6 


47 ± 


40 


poss. broad base 






C III] 


1909 


4644 ± 5 


89 ± 19 


1.2e-18 


9 


220 it 


35 


abs 






C II] 


2326 


Z(\x.Q _|_ 1 9 
jujy zn j-z 


45 ± 10 


3.6e-19 


5 


71 ± 


28 








Mg II 


2799 


6813 ± 17 


146 ± 25 


1.9c-18 


15 


475 ± 110 


abs and atm abs 






[Ne V] 


3426 


8348 ± 5 


50 ± 25 


1.8e-19 


4 


81 it 


46 




5C6.95 


2.877 


Ly a 


1216 


4726 ± 13 


109 ± 9 


5.6c-18 


14 


375 ± 


90 


abs 


Rn9114Q 9-U9Q^fi4fl 
ivy -Duz_i_±4:y ,zxzyou4:U 


-L 


N V 


1240 


A Q1 9 4- Q 


91 4- A 


7 ip in 


1 8 


60 ± 


12 






.007 


C IV 


1549 


5986 ± 20 


92 ± 11 


2.4e-18 


14 


95 ± 


35 


abs 






C III] 


1909 


7396 ± 25 


128 ± 24 


2.7c-18 


14 


117 ± 


40 


abs 


5C6.160 


1.624 


C IV 


1549 


4065 ± 5 


66 it 8 


3.2c-17 


13 


291 ± 


45 


abs 




4- 


He II 


J-O^U 


A'iCiA -t- ^ 


-)- 1 9 

jo m -Lz 


9 7 P 1Q 


1 7 


31 ± 


13 






.002 


C III] 


1909 


5004 ± 7 


78 ± 20 


l.le-17 


6 


165 ± 


45 


abs 






Mg II 


2799 


7351 ± 20 


140 ± 60 


8.0e-18 


1.5 


201 ± 


50 


abs 


5C6.237 


1.620 


C IV 


1549 


4058 ± 15 


70 it 20 


2.6c-17 


12 


78 it 


28 


abs & cr in line 


7C B021749.9+322724 


± 


He II 


1640 


4280 ± 20 


86 it 25 


1.0e-17 


4 


35 ± 


7 






.003 


C III] 


1909 


5007 ± 15 


102 ± 28 


1.4e-17 


5 


52 ± 


12 


abs 






Mg II 


2799 


7332 ± 8 


82 ± 24 


4.0c-18 


2.6 


27 ± 


10 


abs & cr in line 


5C6.251 


1.665 


Si IV / O IV 


1400 


3726 ± 8 


44 ± 8 


5.8c-18 


5 


47 ± 


18 




7C B021850.4+340239 


± 


C IV 


1549 


4107 ± 12 


91 ± 10 


1.2e-17 


8 


108 ± 


15 


abs 




.002 


He II 


1640 


4376 ± 25 


110 ± 60 


3.2e-18 


1.2 


30 ± 


15 








C TTTl 


-Lyuy 


^nfin -t- 90 
ouuu in zu 


114 + 40 


7 9 P 1Q 




85 ± 


24 








C II] 


2326 


6199 ± 6 


72 ± 28 


1.0c-18 


1.3 


21 ± 


9 








[Nc IV] 


2424 


6463 ± 4 


20 it 4 


2.1c-19 


1.1 


4 ± 


3 








Mg II 


2799 


7410 ± 40 


155 ± 25 


5.2e-18 


4 


131 it 


50 


atm abs 


5C6.264 


0.832 


[Oil] 


3727 


6829 ± 2 


23 ± 2 


1.8c-19 


8 


62 ± 


15 




7P Rfl91QzlQ 1 4-llzllzll 


m 






















.001 




















5C6.282 


2.195 


Ly a 


1216 


3902 ± 10 


17 it 5 


3.2e-19 


2.9 


62 ± 


31 


abs 


7C B022140.8+340609 


± 


N V 


1240 


3965 ± 5 


43 ± 12 


4.2e-19 


4 


85 ± 


32 






.003 


C IV 


1549 


4950 ± 20 


150 ± 40 


3.5e-19 


3.5 


71 ± 


32 


abs 






He II 


1640 


5235 ± 4 


80 ± 32 


1.7c-19 


2.9 


27 ± 


16 


abs 






C III] 


1909 


6063 ± 30 


123 ± 40 


1.4e-19 


1.8 


31 it 


11 




5C6.286 


1.339 


C IV 


1549 


3625 ± 7 


74 ± 6 


2.9c-17 


7 


400 ± 100 


abs 


7C B022209.4+314550 


± 


C III] 


1909 


4451 ± 10 


86 ± 15 


6.3e-18 


2.6 


121 ± 


.38 






.003 


Mg II 


2799 


6563 ± 20 


169 ± 25 


8.8c-18 


3.3 


460 ± 100 


abs 
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name 


z 


line 


\ 

Arcst 


Aobs 




FWHM 


nux 


snr 


Wa 




notes 








( A ) 


(A) 




(A) 




m ) 


CT 








5C6.287 


2.296 


Ly a 1 N V 


1216 


4005 


± 


5 


88 


± 


12 


8.5e-17 


15 


340 ± : 


LOO 




7C B022220.0+313901 


± 


C IV 


1549 


5114 


± 


15 


101 


± 


18 


3.5c-17 


8 


192 ± 


80 


abs 




.002 


He II 


1640 


5408 


± 


8 


35 


± 


30 


8.5e-18 


1.4 


57 ± 


35 








C III] 


1909 




1 

± 


12 


oU 


1 

± 


40 


3.7e-18 


1.9 


145 ± 


25 




5C6.288 


2.982 


Ly a 


1216 


4842 


± 


5 


23 


± 


3 


6.1e-19 


12 


240 ± : 


LOO 




7C B022219. 3+310527 


± 


C IV 


1549 


6181 


± 


6 


22 


± 


5 


2.7c-20 


1.4 


36 ± 


15 






.003 


He II 


1640 


6541 


± 


6 


24 


± 


6 


4.2e-20 


1.6 


34 ± 


12 




5C6.291 


2.910 


Lytt/NV 


1216 


4752 


± 


1 


22 


± 


2 


2.9e-18 


90 


372 ± 


60 


abs 


7C B022309.7+340801 


± 


C IV 


1549 


6022 


± 


2 


25 


± 


5 


1.7e-19 


12 


14 ± 


8 


abs 




.002 


He II 


1640 


6417 


± 


2 


30 


± 


4 


2.2e-19 


9 


18 ± 


7 








C III] 


1909 


7465 


± 


2 


145 


± 


25 


7.6c-19 


17 


54 ± 


11 


narrow core 


7C0808+2854 


1.883 


Ly a 


1216 


3518 




12 


105 




20 


5.8c-17 


50 


190 ± 


45 


abs 


7C B080832. 1+285402 


± 


Si IV / O IV 1400 


4020 


± 


10 


117 


± 


25 


5.2e-18 


6 


22 ± 


9 


abs 




.003 


O IV 


lo4y 


4466 


± 


9 


108 


± 


20 


1.9e-17 


15 


90 ± 


18 


abs 






1 1 i l 

Kj 111J 


1909 


5485 


± 


10 


130 


± 


25 


7.0c-18 


8 


46 ± 


11 








Mg II 


2799 


8070 


± 


4 


200 


± 


40 


2.1c-18 


3.3 


24 ± 


11 


abs 


5C7.17 


0.936 


TTTl 
L< 111J 


iyuy 


3695 


± 


15 


95 


± 


25 


1.8c-18 


9 


302 ± 


L40 


abs 


7C B080956.4+270058 


± 


1\/T«- TT 
IVlg 11 


a ( yy 


5407 


± 


15 


202 


± 


60 


1.4e-18 


8 


260 ± : 


L20 






.001 


[in e v j 




6479 


± 


2 


14 


± 


4 


9.6e-20 


1.8 


27 ± 


9 








tti 
[V I1J 


Q707 


7217 


± 


3 


23 


± 


3 


3.5c-19 


9 


76 ± 


18 








flNJo TTTl 

[iNe iiij 


oooy 


7488 


± 


2 


16 


± 


3 


9.9e-20 


2.5 


24 ± 


7 




5C7.70 


2.617 


Ly a 


1216 


4420 


± 


4 


40 


± 


10 


3.9c-19 


8 


88 ± 


37 


abs 


7C B081400.9+292731 


± 


C IV 


1549 


5633 


± 


12 


100 


± 


50 


1.9c-19 


2.6 


37 ± 


25 


abs 




.003 


TJ „ TT 

He 11 


1640 


5947 


± 


14 


140 


± 


70 


1.8e-19 


2.2 


37 ± 


25 








C III] 


1909 


6905 


± 


5 


100 


± 


50 


4.2e-19 


8 


95 ± 


29 


abs 


5C7.85 


0.995 


C III] 


1909 


3799 


± 


7 


70 


± 


30 


9.3e-19 


g 


43 ± 


20 


cr in line 


7C B081445.6+272601 


± 


Mg II 


2799 


5602 


± 


12 


128 


± 


20 


2.1e-18 


10 


77 ± 


22 


abs 




.001 


[O II] 


3727 


7436 


± 


3 


29 


± 


3 


2.6e-19 


3.6 


11 ± 


3 








[Ne III] 


3869 


7718 


± 


5 


29 


± 


3 


2.9e-19 


3.8 


14 ± 


3 




5C7.87 


1.764 


Ly a 


1216 


3370 


± 


4 


28 


± 


16 


3 Oe-18 


2.7 


96 ± 


60 


111 11U15L 


7C B081503.9+245827 


± 


C IV 


1549 


4278 


± 


8 


71 


± 


15 


4.2e-18 


30 


373 ± 


74 


abs 




.002 


He II 


1640 


4534 


± 


5 


69 


± 


23 


4.5e-19 


6 


41 ± 


15 








C III] 


1909 


5263 


± 


10 


52 


± 


30 


9.1c-19 


9 


148 ± 


38 








Mg II 


2799 


7742 


± 


10 


91 


± 


30 


5.8e-19 


3.3 


80 ± 


40 


in atm abs 


5C7.95 


1.203 


C IV 


1549 


3413 


± 


5 


84 


± 


24 


1.7e-17 


21 


226 ± 


57 


r\ b<> on bl n p wi ti o" 


7C B081520. 0+244506 




C III] 


1909 


4195 




8 


78 


-|- 


22 


4.3e-18 


14 


75 ± 


12 






.002 


Mg II 


2799 


6171 


± 


10 


150 


± 


45 


9.5e-18 


16 


194 ± 


46 








[O II] 


3727 


8213 


± 


12 


40 


± 


20 


7.9e-19 


3.1 


22 ± 


10 




5C7.118 


0.527 


Mg II 


2799 


4293 


± 


3 


59 


± 


15 


7.5e-18 


31 


65 ± 


14 


cr in line 


7C B081615. 1+265130 


± 


[Ne V] 


3346 


5109 


± 


2 


19 


± 


4 


8.8e-20 


1.4 


2 ± 


1 






.002 


[Ne V] 


3426 


5230 


± 


2 


15 


± 


3 


9.0e-19 


14 


14 ± 


5 


Ca TT 393*3 abs ? 






H 7 


4340 


6638 


± 


12 


115 


± 


20 


1.2e-18 


5 


29 ± 


f 7 








H/3 


4861 


7424 


± 


8 


120 


± 


18 


3.4e-18 


9 


89 ± 


17 








[O III] 


4959 


7568 


± 


3 


27 


± 


12 


2.1c-19 


2.5 


5 ± 


2 








[O III] 


5007 


7638 


± 


2 


17 


± 


3 


6.4e-19 


8 


17 ± 


5 




5C7.194 


1.738 


Ly a 1 NV 


1216 


3352 


± 


14 


92 


± 


30 


8.9c-17 


32 


183 ± 


37 


abs 


7C B081914.4+254809 


± 


Si IV / O IV 1400 


3841 


± 


6 


97 


± 


18 


l.Oc-17 


7 


37 ± 


11 






.005 


C IV 


1549 


4242 


± 


6 


85 


± 


16 


1.5e-17 


19 


66 ± 


12 


abs 






He II 


1640 


4464 


± 


8 


99 


± 


19 


3.3e-18 


5 


16 ± 


7 








C III] 


1909 


5217 


± 


8 


100 


± 


16 


1.0c-17 


16 


67 ± 


13 








[Ne IV] 


2424 


6637 


± 


5 


35 


± 


5 


5.9c-19 


6 


5 ± 


2 








Mg II 


2799 


7671 


± 


6 


52 


± 


10 


4.0e-18 


25 


42 ± 


10 


atm abs 



Table Al. (cont) 
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name 


z 


line 


Arest 






FWHM 


flux 


snr 


w A 




notes 








(A) 


i A \ 

(A) 




(A) 


(W m~ 2 ) 


a 


(A) 






5C7.195 


2.034 


Ly a 


1216 


3689 ± 


4 


34 ± 4 


3.9c-18 


23 


527 ± 120 




7C B081920. 1+281538 


± 


N V 


1240 


3764 ± 


4 


27 ± 5 


1.4e-18 


11 


183 ± 


70 






.002 


Si IV / O IV 1400 


4241 ± 


15 


40 ± 20 


3.3c-19 


5 


41 ± 


25 


cr in line 






C IV 


1549 


4702 ± 


2 


34 ± 3 


8.5e-19 


17 


95 ± 


40 


abs 






He II 


1640 


4966 ± 


2 


21 ± 4 


3.0e-19 


10 


33 ± 


16 








C III] 


1909 


5772 ± 


5 


41 ± 7 


9.2e-19 


13 


167 ± 


46 


abs 


5C7.230 


1.242 


C IV 


1549 


3468 ± 


6 


51 ± 40 


1.2c-17 


15 


159 ± 


80 


in noise 


7C B082134.3+244829 


± 


He II 


1640 


3674 ± 


5 


52 ± 40 


1.7c-18 


2.8 


26 ± 


6 


in noise 




.002 


C III] 


1909 


4274 ± 


8 


69 ± 17 


2.0e-18 


7 


39 ± 


16 








C II] 


2326 


5222 ± 


2 


25 ± 5 


3.1e-19 


2.6 


7 ± 


3 








Mg II 


2799 


6276 ± 


4 


52 ± 30 


1.5e-18 


7 


40 ± 


9 


abs & cr in line 


7C0825+2930 


2.315 


Ly a 


1216 


4053 ± 


8 


127 ± 18 


9.0c-18 


30 


654 ± 140 


abs 


7C B082505.4+293018 


± 


Si IV / O IV 1400 


4640 ± 


6 


80 ± 22 


9.1c-19 


5 


67 ± 


18 






.001 


C IV 


1549 


5137 ± 


4 


144 ± 24 


5.7c-18 


43 


440 ±110 








C III] 


1909 


6328 ± 


6 


186 ± 32 


2.9e-18 


22 


287 ± 


70 


abs 



Table Al. (cont.) Table summarising the emission lines present in our optical spectra of the 7C quasars and BLRGs. Emission lines 
labelled 'abs' have absorption in them. Emission lines labelled 'atm abs' and 'cr in line' suffer from atmospheric absorption and cosmic 
ray hits, respectively. The line fluxes from the non-spcctrophotomctric observations in Table 2 have been scaled by up to 2.5 mag. to 
account for absorption by cloud. 
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Figure Al. Optical spectrum of the quasar 5C6.264. This object was observed through cloud and the spectrum suffered several 
magnitudes of extinction. The spectrum has a blue continuum and a single narrow emission line which we identify as [Oil] A3727 at a 
rcdshift of 0.832. The if -band image of this object shows a very strong unresolved component and is ~ 1 mag brighter than the mean 
K — z relation for radio galaxies (Eales et al. 1997). An infrared spectrum of 5C6.264 shows a broad Ha line (Paper VIII). We attribute 
the lack of broad lines in the optical spectrum presented here (Mgll A2799 expected at 5130A) as solely due to its poor snr. 
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Figure A2. Optical spectrum of the probable quasar 5C6.288. This spectrum shows strong narrow Lyman-a emission and weak narrow 
CIV and Hell lines. The CIV and Hell lines are clearly visible on the 2-dimcnsional spectrum. The CIII] A1909 line is not observed as it 
falls in the atmospheric absorption band at ~ 7600A. Regions of the spectrum affected by night sky lines and cosmic rays are marked NS 
and CR, respectively. The bright absolute magnitude (Mg = —24.4) and compact nature at ii"-band are suggestive of a lightly reddened 
quasar (Ay ~ 1), which explains the lack of broad UV lines in the optical spectrum (see Paper VIII). 
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Figure A3, (cont) 
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Figure A3, (cont) 
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Figure A3, (cont) Optical spectra of quasars and broad-linc radio galaxies in the 7C sample. 
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APPENDIX B: 3CRR QUASARS 

Table Bl lists the 3CRR quasar sample used in this paper 
which was combined with the 7C sample to give good cover- 
age of the radio luminosity-redshift plane. Only two quasars 
(3C345 and 3C454.3) are excluded on the basis of hav- 
ing Doppler boosted cores/jets which raise their total low- 
frequency flux above the sample limit. One quasar (3C286) 
has a rest-frame 1 GHz spectral index Q ra( j < 0.5, but has 
enough extended flux to remain in the sample. Therefore it 
is called a core-jets quasar (CJS; see Section 2). 

The full list of objects we classify as broad-line ra- 
dio galaxies in the 3CRR sample is 3C22 (RLSE), 3C33.1 
(LJWU), 3C41 (SRL), 3C61.1 (EH), 3C67 (LJWU), 3C219 
(LJWU), 3C268.3 (LJWU), 3C303 (EH), 3C318 (SS), 3C381 
(GO), 3C382 (EH), 3C390.3 (EH) & 3C455 (SDMA). Key 
to broad-emission line references: EH - Eracleous & Halpern 
(1994); GO - Grandi & Osterbrock (1978); LJWU - Laing et 
al. (1994); RLSE - Rawlings et al. (1995); SDMA - Spinrad 
et al. (1985); SRL - Simpson, Rawlings & Lacy (in prep.); 
SS - Spinrad & Smith (1976) 
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Name 


•3178 

(Jy) 


a rad 


z 


lo SlO 1-^151} 
riz sr ) 


V 


Mb 


3C9 


19.492 


1.08 


2.012 


28.68 


18.21 


—26.60 


3C14 


11.336 


0.97 


1.469 


28.10 


20.00 


—24.42 


3C43 


12.644 


0.76 


1.457 


28.09 


20.00 


—23.70 


3C47 


28.776 


0.98 


0.425 


27.39 


18.10 


—23.57 


3C48 


59.950 


0.71 


0.367 


27.46 


16.20 


—25.27 


3C68.1 


13.952 


0.97 


1.238 


28.03 


19.60 


—25.84 


3C109 


23.544 


0.79 


0.305 


26.96 


17.88 


—23.14 


3C138 


24.198 


0.56 


0.759 


27.65 


17.90 


—25.00 


3C147 


65.945 


0.60 


0.545 


27.80 


16.90 


—25.43 


3C175 


19.184 


1.01 


0.768 


27.77 


16.60 


—26.66 


3C181 


15.805 


0.93 


1.382 


28.16 


18.92 


—25.75 


3C186 


15.369 


1.19 


1.063 


27.94 


17.60 


—26.47 


3C190 


16.350 


0.88 


1.197 


28.06 


20.00 


—23.98 


3C191 


14.170 


1.00 


1.952 


28.52 


18.65 


—25.86 


3C196 


74.338 


0.88 


0.871 


28.38 


17.60 


—26.24 


3C204 


11.445 


1.07 


1.112 


27.94 


18.21 


—26.20 


3C205 


13.734 


1.00 


1.534 


28.23 


17.62 


—27.46 


3C297 


14.824 


0.72 


0.684 


27.48 


18.15 


—24.80 


3C208 


18.312 


1.07 


1.109 


28.06 


17.42 


—26.49 


3C212 


16.459 


0.79 


1.049 


27.92 


19.06 


—24.73 


3C215 


12.426 


1.02 


0.411 


26.98 


18.27 


—23.35 


3C216 


22.018 


0.67 


0.668 


27.65 


18.48 


—24.27 


3C245 


15.696 


0.73 


1.029 


27.88 


17.25 


—26.73 


3C249.1 


11.663 


0.88 


0.311 


26.69 


15.72 


—25.32 


3C254 


21.691 


1.05 


0.734 


27.74 


17.98 


—24.75 


3C263 


16.568 


0.87 


0.656 


27.50 


16.32 


—26.22 


3C268 4 


11.227 


0.90 


1.400 


28.01 


18.42 


—26.72 


3C27f) 1 


14.824 


0.80 


1.519 


28.23 


18.61 


—25.54 


3C275 1 


19.947 


0.82 


0.557 


27.52 


19.00 


—23.26 


3C286 


27.250 


0.37 


0.849 


27.58 


17.25 


—25.94 


3C287 


17.767 


0.52 


1.055 


27.79 


17.67 


—26.78 


3C399 1 


24.743 


0.61 


0.904 


27.85 


16.78 


—26.90 


3C334 


11.881 


0.90 


0.555 


27.25 


16.41 


-25.78 


3C336 


12.535 


0.92 


0.927 


27.72 


17.47 


-26.24 


3C343 


13.516 


0.68 


0.988 


27.60 


20.61 


-23.27 


3C351 


14.933 


0.78 


0.371 


26.91 


15.28 


-26.11 


4C16.49 


11.445 


1.00 


1.296 


28.05 


18.50 


-25.62 


3C380 


64.746 


0.72 


0.691 


28.10 


16.81 


-25.91 


3C432 


11.990 


1.07 


1.805 


28.36 


17.96 


-26.61 


3C454 


12.644 


0.79 


1.757 


28.33 


18.47 


-25.86 



Table Bl. Data for the 3CRR quasar sample. The 178 MHz fluxes are those listed in LRL multiplied by 1.09, to correct for an error 
in the KPW flux scale at 178 MHz (Roger ct al. 1973). Spectral indices have been evaluated at 1 GHz rest- frame. Rest-frame 151 MHz 
luminosities arc calculated from corrected 178 MHz fluxes and spectral fitting as described in Blundcll ct al. (1998a). All V magnitudes 
are taken from LRL and have not been corrected for galactic reddening. Absolute magnitudes were calculated as described in Section 
5.1. L151 and Mg assume Qm = 1> = 0- 
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